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PREFACE 


The  Air  Force  SINCGARS  VHF-FM  radio  is  a  frequency-hopping 
anti-jam  radio  which  utilizes  an  electrically  short  antenna  to 
minimize  aerodynamic  drag  on  airborne  platforms.  The  development  of 
optimally  efficient,  electronically  tuneable  antennas  for  this  radio 
is  of  interest.  Although  the  antenna  groundplane  is 
platform-dependent,  it  is  usually  small  compared  to  an  rf  wavelength. 

A  circular  groundplane  provides  a  standardized  groundplane  geometry 
with  which  to  model  and  evaluate  candidate  antennas.  Accordingly,  a 
VHF  antenna  range  witji  an  8  ft.  diameter  circular  groundplane  has  been 
constructed  at  The  MITRE  Corporation  to  evaluate  candidate  antennas. 

The  electrical  properties  of  a  monopole  element  at  the  center  of  a 
circular  groundplane  of  finite  vadius  are  of  interest  to  this  program 
for  (a)  qualifying  the  antenna  range;  (b)  establishing  antenna 
standards  with  which  to  measure  test  antennas;  and  (c)  modeling 
candidate  antennas.  survey  of  the  literature  revealed  that 

although  this  antenna  has  the  simplest  geometry  of  any  monopole 
antenna,  its  properties  are  neither  well  understood  nor  standardized, 
particularly  for  groundplane  radii  which  are  small  or  comparable  to  a 
wavelength.  Tb«  present  paper  attempts  to  address  this  deficiency. 
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M.  Weiner  wrote  Sections  1-5.  S.  Cruze  contributed  to  Section 
3.5  and  wrote  program  MONOPL  in  Appendix  B3.  C.  Li  contributed  to 
Section  3.4;  edited  program  RICHMDl  in  Appendix  B2  and  program 
AWADALLA  in  Appendix  B5;  and  wrote  program  BARDEEN  In  Appendix  Bl, 
program  RICHMD2  In  Appendix  B2,  and  program  MONOSTOR  in  Appendix  B4. 

W.  Wilson  contributed  to  Sections  2.3  and  3.3.  J.  E,  Kriegel  of 
The  MITRE  Corporation  derived  the  correct  form  of  the  continued 
fractions  given  in  Eqs.  3.5-4  and  3.5-5  and  contributed  to  the 
evaluation  of  the  limits  in  Eqs.  3.3-22 — 3.3-24.  W.  C.  Corrleri 
skillfully  performed  the  measurements  discussed  In  Section  5.  K. 
Pamidi  contributed  to  the  development  of  Eq.  3.3-16. 

The  authors  are  grateful  to  Prof.  Alfred  Leltner  of  Rensselaer 
Polytechnic  Institute  for  helpful  conversations  regarding  his  method 
of  oblate  spheroidal  wave  functions;  Prof.  Jack  Richmond  of  Ohio  State 
University  for  helpful  conversations  and  a  magnetic  tape  of  his  method 
of  moments  program;  and  Prof.  Kamal  Awadalla  of  Menoufla  University 
(Egypt)  for  helpful  correspondence  Including  a  listing  of  his  program 
for  the  method  of  moments  combined  with  the  geometric  theory  of 
diffraction. 

After  the  authors  had  obtained  results  using  the  integral  equation 
method  and  the  method  of  oblate  spheroidal  functions,  It  was  possible 
to  confirm  the  correctness  of  Prof.  Richmond's  method  of  moment 
results  which  were  subsequently  published  In  the  open  literature  (see 
Ref.  [2]). 
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SECTION  1 


INTRODUCTION 


Monopole  antennas  are  commonly  employed  In  airborne  and  ground-based 
communication  systems  at  a  wide  range  of  frequencies.  The  electrical 
properties  of  such  antennas  are  dependent  upon  the  geometry  of  both  the 
monopole  element  and  the  groundplane.  Typically,  the  monopole  element  may 
be  electrlcally-short  (length  much  less  than  a  quarterwave)  or 
near-resonant  (length  approximately  a  quarterwave)  and  may  be  thin  (length 
to  radius  ratio  much  greater  than  10^)  or  relatively  thick  (length  to 
radius  ratio  of  10^  to  10^).  In  addition,  the  groundplane  dimensions  may 
vary  from  a  fraction  of  a  wavelength  to  many  wavelengths.  Therefore  it  is 
desirable  to  know  how  the  Input  Impedance  and  radiation  pattern  of  the 
antenna  change  as  the  dimensions  of  the  monopole  element  and  the 
groundplane  vary.  The  directional  gain  on  or  near  the  radio  horizon  (the 
groundplane  is  assumed  to  be  horizontal)  is  of  particular  Interest  since 
the  maximum  operational  range  of  a  communications  system  often  depends  upon 
the  directivity  on  the  radio  horizon. 

This  study  is  restricted  to  a  monopole  geometry  consisting  of  a 
vertical  cylindrical  element  at  the  center  of  a  perfectly  conducting^ 
infinitely  thin,  circular  groundplane  in  free  space.  This  geometry  is  of 
Interest  because  Its  radiation  pattern  Is  uniform  in  the  azimuthal 
direction  and  because  its  electrical  characteristics  are  a  function  of 
primarily  only  three  parameters,  namely,  the  element  length,  the  element 
radius,  and  the  groundplane  radius,  when  each  is  normalized  to  the 
excitation  wavelength.  For  these  reasons  this  geometry  is  conducive  to 
analysis,  experimental  verification,  and  standardization. 


1 


A  typical  feed  for  the  monopole  antenna  Is  a  coaxial  line  whose  Inner 
conductor  Is  connected  through  a  hole  In  the  groundplane  to  the  vertical 
monopole  element  and  whose  outer  conductor  Is  connected  by  means  of  a 
flange  to  the  groundplane.  Typically,  the  Inner  conductor  diameter  Is 
equal  to  the  monopole  element  diameter  and  the  outer  conductor  diameter  Is 
equal  to  the  groundplane  hole  diameter.  Unless  stated  otherwise,  such  a 
feed  will  be  assumed  In  this  study.  The  ratio  of  the  coaxial  line's  outer 
to  Inner  conductor  diameters  affects  the  antenna's  Input  Impedance,  but  only 
significantly  for  a  relatively  thick  monopole  element  on  a  very  small 
groundplane. 

For  the  Idealized  case  of  a  groundplane  of  Infinite  extent  and  of 
Infinite  conductivity,  the  monopole  antenna  may  be  modeled  by  the  method  of 
Images  as  a  dipole  with  one-half  the  Input  Impedance  and  double  the  peak 
directivity  of  the  dipole.  The  Infinite  groundplane  prevents  monopole 
radiation  Into  the  hemisphere  below  the  groundplane  but  allows  a  radiation 
pattern  Identical  to  that  of  the  dipole  In  the  upper  hemisphere.  However, 
for  a  monopole  element  mounted  on  a  groundplane  of  finite  extent,  the  outer 
edge  of  the  groundplane  diffracts  Incident  radiation  In  all  directions  and 
consequently  modifies  the  currents  on  the  groundplane  and  vertical  element 
from  those  of  an  Infinite  groundplane.  At  the  outer  edge  of  the 
groundplane,  the  currents  on  the  top  and  bottom  faces  of  the  groundplane 
are  equal  In  magnitude  but  opposite  In  direction  because  the  net  current 
must  be  zero  at  the  edge.  Outer  edge  diffraction  becomes  Increasingly 
significant  with  decreasing  size  of  the  groundplane  because  of  Increasing 
magnitude  of  the  currents  on  the  groundplane  faces  at  the  outer  edge.  Edge 
diffraction  can  alter  the  Input  Impedance  by  more  than  100%  and  directive 
gain  In  the  plane  of  the  groundplane  by  more  than  6  dB  from  the  values  for 
a  groundplane  of  Infinite  extent. 


Theoretical  models  exist  for  predicting  the  effects  of  diffraction  by 
the  outer  edge  of  the  groundplane.  The  existing  models  may  be  classified 
into  two  categories  distinguished  by  whether  the  current  distribution  on 
the  monopole  element  is  initially  known  or  is  unknown. 

When  the  monopole  element  is  very  thin  and  not  too  long,  its  current 

distribution  is  approximately  sinusoidal  and  independent  of  the  radius  of 

the  groundplane.  Consequently,  the  element's  current  distribution  can  be 

initially  specified  and  only  the  groundplane' s  current  distribution  need  be 

determined.  For  this  category  of  monopoles,  the  theoretical  models 

reported  in  the  literature  consist  essentially  of  Bardeen's  integral 

equation  method  for  the  groundplane  radius  small  compared  to  a 

wavelength^^\  Richmond's  method  of  moments  (groundplane  only)  for  the 

(2) 

groundplane  radius  not  too  large  compared  to  a  wavelength'-  ■' ,  Leltner  and 
Spence's  method  of  oblate  spheroidal  wave  functions  for  the  groundplane 
radius  comparable  to  a  wavelength^^^~^^\  Tang's  scalar  theory  of 
dif f ractlon^^^  and  the  geometric  theory  of  diffraction  (GTD)  for  the 
groundplane  radius  large  compared  to  a  wavelength,  and  Storer's  variational 
method  for  the  groundplane  radius  very  large  compared  to  a 
wavelength^ »(8) ^ 

When  the  monopole  element  is  relatively  thick,  its  current  distribution 
is  no  longer  sinusoidal  and  consequently  the  current  distribution  on  both 
the  monopole  element  and  the  groundplane  need  to  be  determined  as  a 
function  of  the  groundplane  radius.  For  this  category  of  monopoles,  the 
theoretical  models  reported  in  the  literature  consist  essentially  of 
Richmond's  method  of  moments  for  the  groundplane  radius  not  too  large 
compared  to  a  wavelength^^^  and  Awadalla-Maclean's  method  of  moments 
(monopole  element  only)  combined  with  the  geometric  theory  of  diffraction 
for  the  groundplane  radius  large  or  comparable  to  a  wavelength^^^ 

Thiele  and  Newhouse  have  also  reported  a  model  which  combines  the  method  of 
moments  with  the  geometric  theory  of  diffraction  but  their  computer  program 
is  unavailable^^^^ . 


sew 


Each  of  the  existing  models  Is  valid  for  different  and  sometimes 
overlapping  sets  of  limited  values  of  groundplane  radii.  Some  of  the 
models  require  extensive  numerical  computation.  For  these  reasons,  the 
collection  of  models  taken  as  a  whole  has  several  deficiencies.  In  the 
open  literature,  results  for  input  impedance  and  directive  gain  have  never 
been  assembled  as  a  continuous  function  of  groundplane  radius  for  the 
entire  range  of  values  from  zero  to  Infinite  wavenumbers.  In  regions  where 
models  overlap,  it  is  sometimes  unclear  which  models  are  more  accurate.  In 
some  models,  numerical  results  have  been  reported  for  only  a  few  values  of 
groundplane  radius.  In  one  model  (Bardeen's  Integral  equation)  the  base  of 
the  monopole  element  has  not  been  allowed  to  be  in  the  same  plane  as  the 
groundplane  which  is  the  present  case  of  interest.  Computer  programs  are 
not  available  for  some  of  the  older  models  because  the  models  were 
developed  before  the  advent  of  computers.  One  of  the  most  versatile  of  the 
models  (Richmond's  method  of  moments)  gives  only  the  input  impedance  but 
not  the  radiation  pattern.  In  one  model  (Leitner  and  Spence's  oblate 
spheroidal  wave  functions) ,  one  of  the  published  algorithms  for  computing 
the  eigenvalues  is  incorrect.  Finally,  extensive  numerical  results  for 
small  groundplanes  and  for  resonant  monopoles  with  finite  groundplanes  are 
not  found  in  the  open  literature. 

This  paper  attempts  to  correct  these  deficiencies.  Computer  programs 
and  numerical  results  are  presented  for  all  of  the  models.  The  Induced  emf 
method  is  utilized  to  determine  the  input  impedance  of  a  thin  idealized 
monopole  element  in  the  absence  of  a  groundplane.  In  Bardeen's  Integral 
equation  method,  the  excitation  function  for  the  groundplane  currents  is 
extended  to  Include  the  singularity  which  occurs  when  the  base  of  the 
monopole  element  is  in  the  same  plane  as  the  groundplane.  Richmond's 
method  of  moments  is  extended  to  give  the  far-field  radiation  pattern.  In 
Leitner  and  Spence's  method  of  oblate  spheroidal  wave  functions,  the 
published  contlnued'f ractlon  algorithm  for  computing  the  eigenvalues  is 


corrected.  Numerical  results  for  input  Impedance  and  directive  gain  are 
presented  as  a  continuous  function  of  groundplane  radius  for  arbitrary 
radius.  Numerical  results  of  various  models  are  compared  and  the  suspected 
best  available  results  are  Identified.  Extensive  numerical  results  are 
given  for  small  groundplanes  and  for  resonant  monopoles  on  finite 
groundplanes.  New  experimental  data  Is  presented  and  compared  with 
numerical  results. 

Circuit  representations  of  the  monopole  antenna  are  developed  In 
Section  2.  Theoretical  models  and  numerical  results  are  presented  In 
Section  3  for  the  case  In  which  the  current  distribution  on  the  monopole 
element  Is  initially  known.  In  Section  4,  theoretical  models  and  numerical 
results  are  presented  for  the  case  In  which  the  current  distribution  on 
both  the  monopole  element  and  groundplane  are  Initially  unknown.  The 
theoretical  models  are  compared  with  experimental  data  In  Section  5. 
Computer  printouts  of  directive  gain  and  the  far-fleld  elevation  pattern 
are  given  In  Appendix  A.  Computer  programs  of  the  theoretical  models  are 
given  In  Appendix  B. 


SECTION  2 


CIRCUIT  REPRESENTATION 


2.1  Geometry  and  Coordinate  Systems 


Consider  a  monopole  element  of  length  h  and  radius  b  which  is  located 
In  free  space  at  the  center  of  an  Infinitely  thin  circular  groundplane  of 
radius  a  and  of  Infinite  conductivity  (see  figure  1).  The  groundplane 
radius,  when  expressed  in  wavenumbers.  Is  given  by 


£  =  ka 


(2.1-1) 


where 

k  =  271/  X  =•  wavenumber  (m  ^ ) 

X  *  excitation  wavelength  (m) 

The  monopole  element  and  groundplane  have  current  distributions  in  real 
time  given  by 


I(z,t)  -  Re[I(z)eJ‘^‘'],  I(p,t)  -  ReII(P)eJ‘^*') 


(2.1-2) 


where 

u>  =  radian  frequency  of  the  excitation  ••  2  77  c/ X  (rad/sec) 
c  =  wave  velocity  in  free  space  ■  2.9979  x  10  m/s 
I(z),I(P)  “  element  and  groundplane  current  amplitude  distributions, 
respectively  (amp) 

A  field  point  P(r,  0,  0),  expressed  In  spherical  coordinates,  is  shown 
in  figure  1.  The  field  is  uniform  in  the  azimuthal  direction  0.  The 
relationships  between  spherical,  cylindrical  ,  and  oblate  spheroidal 
coordinate  systems  are  shown  in  table  1.  In  the  far-fleld,  the  elevation 
angle  9  Is  related  to  the  oblate  spherlodal  angular  coordinate  T]  by 
77 — ►cos  0  as  the  spherlodal  radial  coordinate  ^  — ►oo. 


TABLE  1 


Spherical,  Cylindrical,  and  Oblate  Spheroidal  Coordinates 


SPHERICAL 

(r,e,(l» 

CYLINDRICAL 
(P.  z,  <» 

OBLATE  SPHEROIDAL 

radius 

r 

P  ^  t  sin  0 

P=*  a[(l-772)(  1+^2)]  1/2 

elevation 

e 

z  *  r  cos  9 

z  =■  a  7?  ^ 

azimuth 

<D 

Note  1:  In  Table  1,  the  notation  ^ Is  that  of  Leitner  and  Spence  (L+S), 
Franklin  Institute  Journal,  Vol.  249,  No.  4,  pp.  299-321,  April  1950.  This 
notation  Is  related  to  that  of  Abramowitz  and  Stegun  (A+S),  "Handbook  of 
Mathematical  Functions",  National  Bureau  of  Standards,  Applied  Mathematical 
Series,  No.  55,  p.  752,  June  1964  by 


^L+S  ^^A+S  ’ 

T]  -  (1  -7?^ 

L+S  '•  'A+S'' 

Note  2;  The  cylindrical  coordinates  (p,  z,  (jl)  are  related  to  the  rectangular 
coordinates  (x,  y,  z)  by 

X  =  p  cos  ^ 
y  »  p  sin 
z  ••  z 


mmmmm 


2.2  Directive  Gain  and  Input  Impedance 


I 


i 


At  a  far-field  observation  point  P(r,  0,  d)),  the  numeric  directive  gain 
d(9,(l))  of  the  antenna  is  defined  as  the  ratio  of  its  radiated  power  density 
s(6,(^)  to  its  radiated  power  density  averaged  over  all  directions.  The 
radiated  power  density  averaged  over  all  directions  is  equivalent  to  the 
power  density  radiated  by  a  hypothetical  "isotropic”  antenna.  Accordingly, 
the  directive  gain  d(0,(l))  expressed  In  spherical  coordinates  with  the 
origin  at  the  antenna,  is  given  by 


d(0,(|l) 


s(0.(ll) 
.27T  /•  TT 


.sn  f  7T 

(1/470  /  I  s(0,(l))  sin  0  d0  d(|) 
Jn  *'o 


(2.2-3) 


For  antenna  patterns  which  are  uniform  in  the  azimuthal  direction,  such  as  for 
the  antenna  geometry  of  figure  l,Eq.  (2.2-3)  reduces  to 


2  a(0) 


2  8'(J7) 


f  s(0)  sin  0  d0  f  s'(r?)  do 

*0  •'-1 


d'(0) 


(2.2-4) 


where 


s(e)  -  (1/2)  Ih^|2  -  IEq|2 

H.  ,  Eq  »  far-field  magnetic  and  electric  field  intensities, 
respectively 

s'(0) .  d'(0)  ”  radiation  power  density  and  directive  gain,  respectively, 
in  oblate  spheroidal  coordinates 

=*  wave  Impedance  in  free  space 

The  numeric  directive  gain  d(0,(|l)  is  related  to  the  directive  gain  D(0,(ll) 
expressed  In  dBi  by 


D(0,(»  -  10  log^^  d(0,(>)  (dBi) 


(2.2-5) 


The  total  time-averaged  radiated  power  of  antenna  is  given 


27r  7T 


total 


If  s(0,(t)  sine  de  d(> 


0  0 


-TT  - 

27Tv^  J  s(0)  sine  de  -  27Tt^  I  s'iV)  drj  (2.2-6) 


The  antenna  radiation  resistance  R,  referred  to  the  complex  amplitude 
l(z“0)  of  the  antenna  base  current,  is  defined  by 


Substituting  Eqs.  (2.2-6)  and  (2.2-7)  into  Eq.  (2.2-4), 


(2.2-7) 


d(0)  =  s(0)  87rr^/[R|  1(2=0)  1^]  -  s'(7))  87nr^/ [rI  I( z=0)  I  =  d'(T]) 


(2.2-8) 


The  antenna  input  Impedance,  is  given  by 


where 


Zin  -  V(z=0)/l(z-0)  =  +  j 


(2.2-9) 


V(z*0)  =  complex  amplitude  of  the  excitation  voltage  across  the 

aperture  of  the  coaxial  line  feed  to  the  antenna  (volts) 

Rj^^  »  input  resistance  (ohms) 

X.  •  input  reactance  (ohms) 

In 


The  input  resistance  R^^  is  related  to  the  radiation  resistance  R  by 


R,  =  R  +  R  .  . 

in  ohmic 


(2.2-10) 


where  R  .  ,  is  the  ohmic  loss  resistance  of  the  antenna  for  finite 
ohmic 

conductivity  of  either  the  raonopole  element  or  the  groundplane.  In  the 

present  paper,  R  .  ,  =0  because  the  monopole  element  and  groundplane  are 

ohmic 

assumed  to  be  of  infinite  conductivity.  Accordingly, 


■'in-"  ■Wo-O  «.2-U) 

Equation  (2.2-11)  Is  a  statement  that  the  antenna  is  assumed  to  have  an 
efficiency  of  unity. 


-.1  '-  I  U-Jv-V  ^  - 


\  A.  V  I 


2.3  Relationship  Between  Radiation  Resistance  and  Directive  Gain 

on  the  Horizon 

For  a  ve^'Clcal  monopole  element  with  a  finite  groundplane,  the 
far-flel'  radiated  power  density  on  the  radio  horizon,  s(0  =7r/2),  is 
determined  only  by  the  current  distribution  on  the  monopole  element  (and 
not  the  groundplane  current)  since  only  the  vertical  monopole  element  has  a 
component  of  electron  acceleration  which  is  normal  to  the  radio  horizon. 
(This  statement  Is  not  true  for  a  groundplane  of  Infinite  extent  since  then 
a  far-fleld  point  on  the  radio  horizon  Is  on  the  groundplane.)  Identical 
monopole  elements  with  Identical  current  distributions,  but  mounted  on 
groundplanes  of  different  finite  radii,  will  consequently  have  Identical 
far-fleld  radiated  power  densities  on  the  radio  horizon.  Accordingly, 

(0  =  7r/2)  •  “  ^/2);  identical  monopole  element.  Identical 

element  current  distribution,  groundplane  of  finite  extent 

(2.3-1) 

where  the  subscript  €  denotes  the  radius  In  wavenumbers  of  the  groundplane 
of  arbitrary  but  finite  extent  (  (  <  oo )  and  the  subscript  0  denotes  a 
groundplane  of  zero  extent  (  €  =0). 

If  Eq.  (2.2-8)  is  substituted  Into  Eq.  (2.3-1)  and  the  quantity 
2  2 

s(0  •  7r/2)  I  l(z’»0)  I  /8  TT  r  Is  computed,  one  obtains  the  following 
relationship  between  radiation  resistance  R  and  numeric  directive  gain  on 
the  horizon  d(0*  Tr/2): 

d^(0  “7T/2)  R^  =■  dQ(9a7r/2)  Rq  »  constant; 

Identical  monopole  elements.  Identical  element  current 
distributions,  groundplane  of  finite  extent  (2.3-2) 


In  Eqs .  (2.3-1)  and  (2.3-2)  the  condition  of  Identical  element  current 
distributions  for  groundplanes  of  different  radii  Is  generally  not 
satisfied  by  monopole  antennas.  The  element  current  distribution  I(z)  Is 
generally  dependent  upon  the  groundplane  current  1(P)  which  In  turn  Is  a 
function  of  the  groundplane  radius.  However,  for  monopole  elements  which 
are  sufficiently  thin  electrically  and  not  too  long,  the  element  current 
distribution  1(2)  Is  approximately  sinusoidal  and  Independent  of  the 
groundplane  current  I(p)  (see  Section  3.1).  Expressions  for  and 

Rq,  for  elements  with  a  sinusoidal  current  distribution,  are  determined  In 
Section  3.2.  Substitution  of  those  expressions  into  Eq.  (2.3-2),  for  the 
case  of  an  infinitely  thin  monopole  element  (b=0),  yields 

d(0  »7r/2)  R  -  (^/47r)  [1  -  cos(kh)]^/sin^(kh);  b=0, 

sinusoidal  element  current  distribution,  groundplanes  of 
finite  extent  (2.3-3) 

where 

T]  »  wave  impedance  of  free  space  =  376.73  ohms 
h  »  length  of  the  monopole  element 

The  condition  b=0  may  be  removed  from  Eq.  (2.3-3)  without  substantially 
altering  the  result  since  the  radiation  pattern  and  radiation  resistance  of 
electrically  thin  elements,  which  are  not  too  long,  are  weakly  dependent 
upon  the  element  radius  (see  Section  3.2). 

2.4  Characterization  of  Currents 

The  characterization  of  the  currents  on  the  monopole  element, 
groundplane,  and  coaxial  line  feed  Is  Illustrated  In  figure  2.  The 
physical  realization,  circuit  representation,  and  two  circuit  idealizations 
of  the  currents  are  shown  In  figures  2(a),  2(b),  2(c),  and  2(d), 
respectively.  In  Fig.  2(b),  the  coaxial  line  excitation  of  Fig.  2(a)  Is 
replaced  by  equivalent  electric  and  magnetic  currents  on  a  conductor 
completely  enclosing  the  coaxial  line.  In  Fig.  2(c),  the  normalized 
ferrite  attenuation  distance  h^/X 1  Is  Idealized  to  be  zero.  In  Fig. 
2(d),  the  magnetic  frill  Is  assumed  to  be  negligible. 


The  currents  of  Interest  are  the  element  current  I(z)  (positive  in  the 

+  z  direction),  the  return  current  I^(z)  on  the  outside  of  the  coaxial  line 

outer  conductor  (positive  In  the  +  z  direction),  the  current  I,  (p)  on  the 

hot 

bottom  face  of  the  groundplane  (positive  in  the  +  P  direction,  and  the 
current  the  top  face  of  the  groundplane  (positive  in  the 

negative  P  direction).  A  net  groundplane  current  I(P)  (positive  in  the 
positive  P  direction)  Is  defined  as 


UP)  ■  -  I,„p(P)  (2.4-1) 

In  the  physical  realization  of  the  currents  [figure  2(a)],  lossy 
ferrite  torlods  are  mounted  along  the  outside  of  the  coaxial  line  outer 
conductor.  Such  a  procedure  is  commonly  utilized  on  antenna  ranges  to 
reduce  the  radiation  or  pickup  of  unwanted  signals  from  currents  induced  on 
the  outside  of  the  transmitter  or  receiver  coaxial  cables.  For  a 
sufficient  number  of  ferrite  toroids  near  the  termination  end  of  the 
coaxial  line,  the  return  current  I^(z)  is  approximately  zero  at  distance 
greater  than  the  1/e  attenuation  distance  h^  from  the  termination  end  of 
the  line.  In  this  pa^^er,  the  ferrite  toroids  are  assumed  to  be 
sufficiently  lossy  so  that 


h^  «  A  (2.4-2) 

where  A  Is  the  excitation  wavelength. 

The  radii  of  the  outer  and  Inner  conductors  of  the  coaxial  line  are  b^^ 
and  b,  respectively,  where  b  is  also  the  radius  of  the  monopole  element. 

The  wall  thickness  of  the  coaxial  line  outer  conductor  Is  assumed  to  be 
much  less  than  Its  diameter.  Consequently,  the  return  current  I^(z)  occurs 
at  the  radial  coordinate  p  ~  b^ , 


15 


The  constraints  on  the  various  currents  are: 


m 


element 

groundplane 

coaxial  line 

element-groundplane 

groundplane-coaxlal  line 


I(2-h)  -  0 
I(p=a)  »  0 

0,  -oo<  z  <-h^ 
l(z»0)  =  -I(P»b) 


(2.4-3) 

(2.4-4) 

(2.4-5) 

(2.4-6) 

(2.4-7) 


The  element  and  groundplane  constraints  are  a  consequence  of  an  open 
circuit  at  the  end  of  the  element  and  the  groundplane.  The  coaxial  line 
constraint  Is  a  consequence  of  the  lossy  ferrites.  The  element-groundplane 
constraint  Is  a  consequence  of  conservation  of  charge  (Klrchhoff's  cui  rent 
law)  at  a  node.  The  groundplane-coaxlal  line  constraint  Is  a  consequence 
of  conservation  of  charge  along  a  conductor. 

(12) 

By  the  use  of  the  equivalence  principle  the  coaxial  line  feed 
excitation  may  be  replaced  by  equivalent  tangential  field  excitations 
defined  along  a  surface  completely  enclosing  the  coaxial  line.  At  field 
points  external  to  this  surface  the  equivalent  field  excitations  will  give 
the  same  field  as  the  original  source  excitation.  In  the  circuit 
representation  of  the  monopole  antenna  currents  [figure  2(b)],  the  coaxial 
aperture  excitation  Is  replaced  by  an  equivalent  surface  magnetic  current 
density  (magnetic  frill)  which  sits  on  top  of  a  thick  groundplane  of 
radius  b^ .  The  magnetic  frill  Is  defined  to  be  positive  In  the  positive 
azimuthal  direction  and  Is  given  by 


-  V(0)/Ipin(b^/b)],  b<p<bj 


,  p<  b,  P  >  b, 


(2.4-8) 


Eq.  (2.4-8)  is  derived  as  follows.  The  radial  field  E  of  the  coaxial 
line  aperture,  assuming  a  TEM  mode  excitation,  is  given  by 


V(0)/pln(bj^/b),  b<P<h^ 


P<  b 


(2.4-9) 


where  V(0)  is  the  positive  voltage  at  z*0  across  the  aperture  with  the 

coaxial  outer  conductor  at  zero  potential.  By  the  equivalence 

(12) 

principle  ,  an  aperture  field  may  be  replaced  by  a  magnetic  frill  M 
which  sits  on  top  of  a  groundplane  congruent  with  the  aperture  surface  and 
which  is  defined  as 


M  =  ,  X  n 

tangential 


(2.4-10) 


where  n  is  the  outward  normal  to  the  aperture  surface  and  E  ,  ,  is 

tangential 

the  tangential  field  at  the  aperture  surface.  Substituting  Eq.  (2.4-9) 
into  Eq.  (2.4-10), 


M  -  E^)  X  u^ 


(2.4-11) 


where  Mx  is  given  by  Eq.  (2.4-8).  Eq.  (2.4-8)  agrees  with  the  result 

(2) 

obtained  by  Richmond  . 

In  the  circuit  representation  of  figure  2(b),  the  net  groundplane 
current  I(p)  is  the  same  as  defined,  by  Eq.  (2.4-1)  with  the  additional 
current  constraint 


0<p<b^ 


(2.4-12) 


i 


If  the  circuit  representation  of  the  monopole  antenna  in  figure  2  is 
now  idealized  by  setting  the  ferrite  1/e  attenuation  distance  hj^/A  =  0,  then 


I^(z)kO  for  z  <  0,  h^/X=  0 


(2.4-13) 


Consequently,  the  coaxial  line  outer  conductor  may  be  removed  from  the 
circuit  as  shown  in  figure  2(c).  The  groundplane  -  coaxial  line  current 
constraint  of  Eq.  (2.3-7)  is  not  disturbed  by  such  an  idealization.  Since 
it  has  already  been  assumed  that  hj^«X,  the  idealization  hj=  0  does  not 
appreciably  alter  either  the  radiation  pattern  or  the  input  Impedance  of 
the  monopole  antenna  provided  that  the  monopole  length  h  »  h^.  Finite 
currents,  on  an  aperture  which  is  small  compared  to  the  excitation 
wavelength,  contribute  little  to  the  far-field  and  input  Impedance 
(compare  with  the  results  for  an  electrlcally-small  dipole).  Experimental 
radiation  patterns  and  measurements  of  input  Impedance  (see  Section  5) 
confirm  that  the  use  of  lossy  ferrite  toroids  along  the  coaxial  line  outer 
conductor  yields  results  which  are  in  close  agreement  with  theoretical 
results  for  the  circuit  Idealization  condition  of  Eq.  (2.4-13)  —  even  for 
electrlcally-small  groundplanes . 

In  the  idealization  of  the  monopole  antenna  circuit,  the  magnetic  frill 
may  be  removed  [figure  2(d)]  without  appreciably  altering  the  radiation 
pattern  or  input  Impedance,  provided  that 


kb^  «  1 

I(z)  or  I(p)  is  initially  known 


conditions  for  neglecting  magnetic  frill 

(2.4-14) 


' 


The  condition  kb^  <<  l  corresponds  to  the  condition  for  TEM  mode  excitation 
of  the  coaxial  line  and  for  negligible  power  being  radiated  from  the 
coaxial  line  aperture.  If  either  l(z)  or  I(P)  Is  initially  known  then  the 
one  which  is  not  known  may  be  determined  from  the  other  without  requiring  a 
knowledge  of  the  original  coaxial  line  excitation  or  its  equivalent  — 
provided  that  the  field  radiated  by  the  known  current  distribution  is  the 
predominant  field  incident  on  the  conductor  whose  current  distribution  Is 
unknown.  When  neither  1(2)  nor  I(P)  Is  known,  then  the  original  source 
.xcltatlon  or  It,  equivalent  (In  thl.  eaaa,  the  magnetic  frill  be 

specified  to  determined  the  unknown  radiated  fields. 


SECTION  3 


MODELS  IN  WHICH  THE  CURRENT  DISTRIBUTION  ON  THE  MONOPOLE 
ELEMENT  IS  INITIALLY  KNOWN 


3.1  Boundary  Conditions 

The  current  amplitudes,  I(z)  and  I(P),  on  the  monopole  element  and 
groundplane,  respectively,  are  generally  complex  and  Initially  unknown 
quantities  (see  figure  1).  Consider  now  the  case  where  the  current 
distribution  on  the  monopole  element  Is  assumed  to  be  sinusoidal.  For  such 
a  case  and  for  the  waveform  dependence  given  In  Eq.  (2.1-2), 


1(2)  =  71^^^)  sin[k(h-z)],  0  £  z  <  h 


(3.1-1) 


where 

h  =  monopole  length  (m) 
k  =  2  TT  /  X  =  wavenumber  (m  ^) 

1(0)  =  current  amplitude  of  the  monopole  base  current  at  z*0  (amp) 

From  Eq.  (3.1-1) 

arg  I(z)  =  constant,  0  £  z  £  h  (3.1-2) 

Although  a  sinusoidal  distribution  of  current  is  not  possible  even  for 
an  Infinitely  thin  antenna,  Eq.  (3.1-1)  is  most  likely  a  fair  approximation 
to  the  current  If  the  monopole  element  Is  sufficiently  thin  electrically 


and  not  too  long.  '  For  a  center-fed  dipole  of  radius  b  and  total 
length  2h,  Elllot^^^^  gives  examples  where  the  current  distribution  Is 
approximately  sinusoidal  and  Is  of  approximately  constant  argument  for 
b/ X  *  1.0  X  10  ^  and  h/ X  ”  0.125,  0.25.  However  for  b/X  “  1.0  x  10  ^  and 
h/ X  “  0.375,  0.5  Elliot  demonstrates  that  the  current  distribution  Is  no 
longer  sinusoidal  near  the  center  of  the  dipole  nor  Is  arg  z  approximately 
constant.  Balanls^^^^  shows  that  for  b/X=  2.0  x  10  ^  and  h/X  =  0.25,  0.5 
the  current  distribution  Is  not  sinusoidal  near  the  center  of  the  dipole. 
Elliot  and  Balanls  also  demonstrate  that  for  b/X  =  1.0  x  10  and  h/X  = 
0.125,  0.25,  0.375,  0.50  the  current  distribution  Is  neither  sinusoidal  nor 
of  constant  phase  and  that  the  deviations  from  Eqs .  (3.1-1)  and  (3.1-2) 
increase  with  Increasing  values  of  h/X  and  b/X.  On  the  basis  of  the  above 
results.  It  appears  that  Eqs.  (3.1-1)  and  (3.1-2)  are  approximately  valid 
for  the  conditions 

b/X  <  10“^,  h/X  <  0.25  (3.1-3) 

In  addition  to  the  constraint  on  I(z)  given  by  Eq.  (3.1-1),  assume  that 
the  return  current  1^(2)  on  the  outside  of  the  coaxial  line  outer  conductor 
(see  figure  2)  Is  given  by 

lj.(z)«0  ,  z  <  0  (3.1-4) 

The  constraint  of  Eq.  (3.1-4)  corresponds  to  the  Idealized  condition  that 
the  ferrite  toroids  have  a  1/e  current  attenuation  distance  h^  given  by 

h^/X=  0  (3.1-5) 

It  should  be  noted  that  Eqs.  (3.1-4)  and  (3.1-5)  do  not  alter  the 
constraint  l^(z=0)  =  ^bot^”  ^1^  given  by  Eq .  (2.3-7)  nor  do  they  Impose 
any  constraints  on  ^bot^’’  ^1^  where  ^bot^*  ^1^  current  on  the 

bottom  of  the  groundplane  at  a  radius  equal  to  that  of  the  outer  conductor. 


Combining  the  current  constraints  given  by  Eqs.  (3.1-1)  and  (3.1-4) 
with  those  given  by  Eqs.  (2.4-3)  -  (2.4-7)  and  Eq.  (2.4-12),  the  current 
constraints  on  the  monopole  antenna  are  given  by 


element 

groundplane 


coaxial  line 
element-groundplane 


I(z)  =  [I(0)/8ln(kh)]sln[k(h-z)],  0  <  z  <  h  (3.1-6) 


I(p=a)  =  0 

lj.(z)~0,  z  <  0 
l(z=0)  =  -I(P*b) 


groundplane-coaxlal  line  I,  (p)  -  I  (z*0) 
-  bot  r 


(3.1-7) 

(3.1-8) 

(3.1-9) 

(3.1-10) 

(3.1-11) 


Tn  Section  3,  it  will  be  assumed  that  all  models  satisfy  the  current 
constraints  given  by  Eqs.  (3.1—6)  —  (3.1-11).  The  results  are  expected  to 
be  approximately  correct  If  the  monopole  element  Is  electrically 
sufficiently  thin  and  not  too  long  [Conditions  (3.1-3)]  and  if  the  ferrite 
toroids  are  sufficiently  lossy  [hj^/X«  1»  Idealized  by  Condition  (3.1-5)]. 
For  these  conditions,  the  circuit  representation  of  the  monopole  antenna  Is 
shown  In  figure  2(d). 

For  the  current  constraints  of  Eqs.  (3. l-6)-(3 . 1-11 ) ,  the  total 
magnetic  and  electric  field  intensities  ^(total)  arbitrary 

field  point  P(x,y,z)  external  to  the  element  excitation  source  points  are 
simply  the  vector  sum  of  the  fields  resulting  from  the  element  current  and 
the  current  Induced  on  the  groundplane  by  the  fields  Incident  from  the 
element. 

Accordingly, 

^(total)  .^(e)  4.-ff<8)ff<e)^^(e)j  .  current  constraints  of 


-E(total)  .-^(e) 


,  Eqs.  (3.1-6)-(3.1-ll) 


(3.1-12) 


where 


H  ,  E  =  magnetic,  electric  field  Intensities,  respectively,,  generated 
by  the  element  current 

^(g)  ^(g) 

H  ,  E  °  =  magnetic,  electric  field  Intensities,  respectively,  generated 

by  the  groundplane  current  Induced  by  the  element  incident 
fields . 

The  element  fields  H  ,  E  are  determined  In  Section  3.2.  The 

-^(o)  -^(o) 

groundplane  fields  H  ,  E  are  determined  In  Sections  3.3  -  3.8  for 
groundplane  radii  of  various  extents. 

3.2  Induced  EMF  Method,  Groundplane  of  Zero  Extent 

Concept  of  a  groundplane  of  zero  extent 

Consider  a  monopole  antenna  excited  by  a  coaxial  line  whose  outer 
conductor  of  radius  b^  is  terminated  by  free  space  rather  than  by  a 
groundplane  [figure  3(a)].  The  groundplane  for  such  an  antenna  is  denoted 
as  being  of  "zero  extent".  As  was  shown  In  Section  2.4,  the  coaxial  line 
excitation  may  be  replaced  by  an  equivalent  magnetic  current  (frill) 
sitting  on  top  of  a  thick  groundplane  of  radius  P=  b^^  [figure  3(b)].  For 
sufficiently  lossy  ferrite  toroids  along  the  outside  of  the  coaxial  line, 
the  current  on  the  exterior  of  the  coaxial  line  outer  conductor  may  be 
neglected  [figure  3(c)].  The  magnetic  frill  may  be  removed  from  the 
circuit  without  appreciably  affecting  the  results  since  kb^  <<  1  for  the 
assumed  sinusoidal  current  distribution  on  the  monopole  element  [see  Eq. 
(3.1-3)  and  the  discussion  concerning  the  circuit  Idealization  of  figure 
2(d)].  Finally,  the  groundplane  of  radius  b^  may  be  removed  from  the 
circuit  without  appreciably  affecting  the  results  since  a  finite  current  on 
an  electrlcally-small  conductor  does  not  radiate  appreciable  power  compared 
to  the  power  radiated  by  a  monopole  element  of  length  h  »  b.^^.  The  circuit 
Idealization  of  a  monopole  antenna  with  a  groundplane  of  zero  extent  Is 
therefore  an  electrically  thin  monopole  element  with  no  groundplane  [figure 


Characterization  of  Currents  for  Monopole  Antenna  with  Groundplane  of  Zero  Extent 
(a)  Physical  realization 


The  near-fields,  far-flelds,  and  input  impedance  of  an  electrically 
thin  monopole  element  are  derived  and  summarized  in  the  remainder  of 
Section  3.2.  The  input  Impedance  is  derived  by  the  Induced  emf  method. 

Near-fields 


Consider  a  monopole  element  of  length  h  with  a  sinusoidal  current 
distribution  I(z')  =  [l(0)/sln(kh)]  sin  fk(h-z')].  0  £  z'  _<  h,  at  points 
Q(x',  y' ,  z')  on  the  surface  of  the  element  (figure  4).  For  an 
electrlcally-thln  element  with  a  known  current  distribution,  the  fields  at 
arbitrary  points  P(x,  y,  z)  external  to  the  element  may  be  determined 
almost  exactly  be  approximating  the  source  points  to  lie  on  the  element 
axis,  l.e.,  Q(x',  y' ,  z')wQ(0,  0,  z').  For  the  current  waveform  of 
Eq.  (2.1-2),  the  magnetic  vector  potential  A  is  given  by^^^^ 


A(x,y,z)  =*  u  A 
z  z 


/\ 

u 


UMO) 


z  47rsln(kh) 


_ 

/  sinfk(h-z^)1  expf-lkPQ)  dz' 

PQ  (3.2-1) 


where 


-7 


o 

PQ 


permeability  of  free  space  =  47r  x  10  henries/m 

2  2  2  1/2 
[(x-x')^  +  (y-y')^  +  (z-z')^]"'^ 

r  2  .  2  ^  ,  ,.2,1/2  .2  ^  ,  ,.2,1/2 

[x  +  y  +  (z-z  )  ]  =  [p  +  (z-z  )  ] 


^  /N 

Up,  U|j,u^  =  unit  vectors  in  the  radial,  azimuthal,  and  axial 
cylindrical  directions,  respectively. 


The  magnetic  and  electric  field  intensities,  H  and  E,  respectively,  are 
given  by 


(7x  A) 


r  1 

“(Jl 


H  - 


1 

Wo- 


(3.2-2) 


jtje 


(7x  H) 


/V 

•  u  , 


p  iO)€^  ^  Z 


/N 

+  u 


1 


z  ja>€ 


_1_ 

P 


-12 


hP 

(3.2-3) 


(PH.) 


where  e  =  permittivity  of  free  space  =  8.854  x  10  farad/m 


Exact  closed  form  expressions  of  A,  H,  E  and  radiation  resistance  were 
first  obtained  for  an  infinitely  thin  element  by  Brlllouln^^^^  and  are 
summarized  by  Stratton''  .  The  magnetic  and  electric  fields  are  given  by 


H  =  J _ _ 

(tl  47rpsl 


1(0)  r 

teln(kh) 


exp(-jk.r^)  -  cos(kh)  exp(-jkr) 


- jz  sln(kh)  exp(-jkr) 

r 


(3.2-4) 


J  H0)r] 


(z-h)  exp(-jkr^) 


4npsin(kh) 


z  cos(kh)  exp(-jkr) 
r 


^  sln(kh)  ^  z  exp(-jkr) 

^z  kr 


J  1(0) n  (z-h)  exp(-jkr^)  ^  ^  co8(kh)  exp(-jkr) 
4npsin(kh)  r  ”  r 


sln(kh)  expC-jkr)!*  1  z^  _  jz^"! 

’  kr'^  '  r^J 


(3.2-5) 


~j  1(0)^  I  ^^P(  ^  cos(kh)  exp(- jkr) 


4  TT  sln(kh) 


,  sln(kh)  B 

az 


exp(- 

r 


o 

-jkr) 


-J  I(0)?7 


4;rsin(kh)  j 
z  sln(kh)  exp(-jk 


exp(- jkr  )  .  X  /.IX 

o  _  cos(kh)  exp(- jkr) 


(3.2-6) 


28 


■  wave  Impedance  of  free  space 

IP^  + 

(P^  +  z2)l/2 


Eqs .  (3.2-4)  -  (3.2-6)  are  Identical  to  the  results  given  by  Stratton 
[^2  ”  ”  h,«=  kh,  r2  =  r^,  r^  =  r,  =  -I(0)/(8ln  (kh)] 

after  the  substitution  of  -J  for  j  to  account  for  the  exp  (-jo;t)  waveform 
of  Stratton  instead  of  the  exp(jaJt)  waveform  of  Eq .  (2.1-2).  The  fields 
given  by  Eqs.  (3.2-4)  -  (3.2-6)  are  exact  for  an  inf Inltely-thin  element 
and  are  almost  exact  for  an  electrlcally-thln  element  with  the  same 
sinusoidal  current  distribution. 

Far-f ields 

Consider  a  field  point  P  at  a  sufficiently  large  radial  distance  r 
which  satisfies  both  the  far  radiation  zone  and  Fraunhofer  diffraction 
conditions  given  by 


h  «  r,  kr  >>  1;  far  radiation  zone  conditions 
kh  /2r  \<  271;  Fraunhofer  diffraction  condition 


(3.2-7) 

(3.2-8) 


For  these  conditions, 


exp(- jkr^) S3  exp[-kj( r-h  cos0)],  h  «  r,  (kh‘^/2r)  «  2  77 


1/r  ss  1/r,  h  «  r 
o 


z  -  haz  =  r  cos0,  h  <<  r 


(1/kr)  -  (z^/kr^)  -  (jz^/r^)ss-j  cos^0,  kr  »  1 


(j/r^)  +  (l/kr^)a  (j/r^),  kr  »  1 


(3.2-9) 


For  the  approximations  of  Eq.  (3.2-9),  the  "exact"  flelrls  given  by 
Eqs.  (3.2-4)  -  (3.2-6)  reduce  to  the  far-flelds  given  by 


H  =  j  T(0)  exp(-ikr) 

(!)  47Tr  slnft  sln(kh) 


exp  ( ikhcosA)  -  cos(kh)  -  j  cosP*  sln(kh) 

(3. ->-10)  -I 


i  I(n)^cosfl  exp(-ikr) 
47rr  slnft  sln(kh) 


jl  -  cos(kh)  -  i  cosP*  sln(kh)J 

(3.2-11) 


>] 


.i  1(0)  ^  exp(-  ikr) 


4  TT  r  sln(kh) 


1  -  cos(kh)  -  i  cosft  sln(kh) 

*-  -1(3. '1-1'’) 


The  resultant  electric  field  F  =u  F  +  u  E  reduces  In  the  far  field  to 

P  p  z  z 

r3.2-13) 


E  far-fleld 


where 


,  j  1(0)??  exp(-  ikr) 
0  47rr  sln0  sln(kh) 


1  -  cos(kh)  -  i  COS0  sln(kh) 


The  time-averaged  Poynting  vector  S  =  (1/2)  (F  x  H)  = 

/\ 

(1/2)  (Uq  F^  X  u^  Hj||)  =  u^  s(0)  with  a  time-averaged  radiated  newer  densltv 
s(0)  given  by 


where 


s(0)  =  lF„r/(2l?) 


3177  r  f 


si  fkh) 


(3. 7-1 4) 


2 

(cos(kh  cos<3)  -  cos(kh)l 
r  ( H )  = - - — 

sin  9 


n 

+  rsln(kb  COS0)  -  cos  P  sln(kh)l 


The  direction  of  maximum  radiated  power  Is  9  =  77/2  rad. 


Substituting  Eq.  (3.2-14)  Into  Eqs.  (2.2-6)  and  (2.2-7),  the  radiation 
resistance  R,  referred  to  the  base-current,  Is  given  for  an  Infinitely  thin 
element  by 


r 

n  Jq  f(e)  sine  de 

R  =  - A - (ohms) 

8  IT  8ln^(kh) 


- ^ -  [Cln(2kh)  -  8ln^(kh)],  b=0 

ATTsln  (kh) 


(3.2-15) 


where  Cln(z)  =  modified  cosine  Integral 


(19) 


ii 

t 


(1-cos  t). 


The  result  given  for  the  definite  Integral  of  Eq.  (3.2-15)  Is  readily 
obtained  by  letting  t  =  cos  0  and  noting  that 


where  x  =  1+t  and  y  =  1-t. 


Substituting  Eqs.  (3.2-14)  and  (3.2-15)  Into  Eq.  (2.2-8),  the  numeric 
directive  gain  d(9)  Is  given  by 


d(9)  = 


_ f(e) 

Cln(2kh)  -  sln^(kh) 


(3.2-16) 


where  f(9)  Is  defined  In  Eq.  (3.2-14). 


For  electrlcally-short  monopole  elements  (kh  «  1)  and  376.73 

w  120  77  ohms,  Eqs.  (3.2-14)  -  (3.2-16)  reduce  to 


f(6).  1 

kh  «  1 

(3.2-17) 

8(9)  -  ‘5 

|2  (kh)2  sin'e  J 

(3.2-18) 

16 

rrr 

R  =  5  (kh)^  (ohms),  kh  «  1,  b  =  0 

(3.2-19) 

d(0)  »  ®, 

kh  «  1 

(3.2-20) 

The  relative  power  radiation  pattern  s(0)/s  =  s(e)/s(7r/2) ,  radiation 

lUdX 

resistance  R,  and  the  directive  gain  (directivity)  d(0)  are  tabulated  in 
Appendix  A1 ,  Tables  Al-1 ,  Al-14,  Al-26,  Al-38  for  h/X  =  1/4,  1/10,  1/40, 
and  1/100,  respectively.  The  numeric  directive  gain  pattern  d(0)  for  a 
quarterwave  monopole  element  [see  figure  8(a)  in  Section  3.9]  Is  similar  to 
that  of  a  half-wave  dipole  except  that  its  peak  directivity  is  less  (1.88 
dbl  vs.  1.76  dBi)  and  its  3  dB  beamwidth  is  more  (94  deg  vs.  78  deg). 

Input  Impedance 


The  input  impedance  may  be  determined  by 
Introduced  by  Brillouln^^^\  With  reference 
impedance  is  given  by^^^^ 


the  Induced  emf  method 
to  figure  4,  the  input 


Z 

in 


b. 


z 


where 


z')  E  (P  =  b,  z  =  z')  dz' ,  b  >  0 
z 

(3.2-21) 


I  =  l(0)/sln(kh) . 
m 


32 


'.VYj 


The  Induced  emf  method  Is  Indeterminate  for  collinear  current  elements 
unless  the  elements  are  of  radius  b  >  0.  For  sufficiently  thin  elements, 
the  current  is  given  by  Eq.  (3.1-1)  and  the  electric  near  field  is 
given  almost  exactly  by  Eq.  (3.2-6).  Substituting  Eqs.  (3.1-1)  and  (3.2-6) 
into  Eq.  (3 . 2-21 ) , 


'in  4irsln(kh) 


dz'  sln[k(h-z')]i 


exp(-jkr  ) 


cos(kh)  exp(-jkr)  .  sin(kh)  f  exp(-jkr)  1 

r  ^z'  L  r  J 


(3.2-22) 


where 


-  u^2,l/2 

r  =  [b  +(z  -h)  ] 

°  ,^2  ^  .2,1/2 

r  “  (b  +  z  ) 


Each  of  the  three  terms  of  the  integrand  of  Eq.  (3.2-22)  may  be  integrated 

/  lO  \ 

by  the  methods  summarized  by  Stratton'  .  However,  the  third  term  of  the 
integrand  should  first  be  integrated  by  parts. 

Accordingly, 


Z,  =  R,  +  J  X,  =  Z,  +  Z,  +  Z, 
in  in  in  1  2  3 


(3.2-23) 


where 


Airsin  (kh) 


J  dz'  sin[k(h-z)]  ^  +  jX 


1  -’“1 

(3.2-2A) 


- ^  Cin(x, )  +  — ^  Cin(x,)  -  Cin(x-) 

Ansin  (kh) 


(3.2-25) 


- \ -  ~  Si(x  )  + -^  Si(x  )  -  Si(x  )  (3.2-26) 

Ansin^(kh)  L  ^ 


j  y?  r  -dz'  sin[k(h-z')]  co8(kh)  exp(- jkr)  *  R2+  ji 

4,,ln^kh)  ^  '  (3.2-27) 


4jTSln  (kh) 


r  8tn(2kh)  _  si(x2)] 


+  cos^(kh)  (-Cln(Xj)  +  Cin(x2)  +  Cln(Xj^)] 


(3.2-28) 


ri 


4nsln  (kh) 


- [cin(x2)  -  Cln(Xj^)  +/n(x^/x2)] 


+  [Sl(x^)  +  Si(x2)  -  2  SKx^)] 


(3.2-29) 


4iTsln  (kh 


n 

rdz'  sin[k(h-z’)]  sln(kh)  exp(-jkr) 

'C\  /I 


(3.2-30) 


sin  (kh) 


4Ttsin  (kh) 


|cin(Xj^)  +  Cin(x2)  -  2  Cin(x2) 


2(Xi  +  x^)  (sin  x^  +  sin  x^) 


4  x^  sin  x^ 


2  2  2 
x^  +  X2  +  2(kb) 


2  2 
X3  +  (kb)^ 


sln(2kh) 


Si(x2)  -  Sl(Xj^)  - 


2(Xi  +  X2)  (cos  x^  -  cos  x^) 


2  2  2 
x^  +  X2  +  2(kb)^ 


(3.2-31) 


n 


4nsln  (kh) 


sin  (kh) 


Sl(x^)  +  Si(x2)  -  2Si(x2) 


2(Xj^  +  X2)  (cos  Xj^  +  cos  X2)  4  x^  cos  x^ 


2  2  2 
x^  +  X2  +  2(kb)^ 


Xj  +  (kb) 


sin(2kh) 


-Cln(Xj^)  +  Cln(x2)  + 


2(Xi  +  X2)  (sin  x^  -  sin  X2) 


2  2  2 
X.  +  X-  +  2(kb) 


(3.2-32) 


-  k[(b  +  h  )  '  +  h] 

X2  -  k[(b^  +  -  h] 

»  kb 


Cln(x)  =  modified  cosine  Integral^^^^  *  f  — (1-cos  t) 

•'o  ^ 

Si(x)  =  sine  Integral^^^^  =  dt 


Summing  Eqs.  (3.2-23)  -  (3.2-32), 

ri  L.  , 


._  *  - 5 -  Cln(x  )  +  Cin(x  )  -  2  Cin(x  ) 

4nsin'^(kh)  ^ 


,  8ln(2kh) 

(Xi  +  X2)(cos  -  cos  X2)  ■ 

2 

2  .  2  ...  2 
x^  +  X2  +  2x^ 

,  [■  (x.  +  x-)(8in  X.  +  sin  x,) 

+  8ln‘^(kh)  — - - - 5 - — - — 

L  «i  +  xj  2*3 


C.  =  - = -  Si(x.)  +  Sl(x,)  -  2  Sl(x,) 

4rrsln2(kh)  ^  ^  3 


sin  X. 


(3.2-33) 


sln(2kh) 

2 


(x^  +  X2)(sln  x^  -  sin  x^) 

2  2  2 

xj  +  x^  +  2x; 


„  r (x,  +  x„)(cos  x^  +  cos  x„)  cos  x» 

+  sln^kh)  - 2 - ^“2 - ^ - ^ 

Xj  ■*■  ^^2  ^  ^*3  ^ 

I 

2  2  1/2  2  2  1/2 
where  x,  -  k[(b  +  h  )  '  +  h),  x„  =  kr(b  +  h  )  '  -  h],  x,  =  kb. 


(3.2-34) 


fl 


For  b=0,  Eqs.  (3.2-33)  and  (3.2-34)  reduce  to 


47Tsln^(kh) 


Cin(2k.h)  -  sin  (kh) 


,  b=0 


in 


Artsin  (kh) 


Si(2kh)  - 


sin  (kh) 
kb 


,  b-.0 


(3.2-35) 

(3.2-36) 


Eq.  (3.2-35)  agrees  with  the  result  for  radiation  resistance  given  by 

Eq.  (3.2-15).  A  comparison  of  Eq.  (3.2-33)  with  Eq.  (3.2-35)  reveals  that 

the  input  resistance  is  relatively  insensitive  to  the  monopole  element 

radius  b  for  kb  <<  1.  The  input  resistance  given  by  Eq.  (3.2-35)  is 

plotted  in  figure  5.  For  kh  =7r/2  radians,  R  =  19.4  ohms. 

in 

The  input  reactance  is  sensitive  to  the  element  radius  b  as  seen  in 
Eq .  (3.2-36).  For  b  «  0,  the  input  reactance  is  given  by 


-  00  ohms ,  kh?*n7r,  n  =  1,  2,  3  ... 
ooohmc,  kh  =  nTT,  n  *=  1,  2,  3  ... 


b  =  0 


(3.2-37) 


From  Eq .  (3.2-36),  resonance  (X  =  0)  occurs  for 

in 

(kh)  s:n7r+  [kb  Sl(2n7r) ] ^ ,  n  =  1,  2,  3  ...  (3.2-38) 

resonance  — 

The  input  resistance  at  resonance,  for  kh  given  by  Eq.  (3.2-38),  is  found 
from  Eq .  (3.2-35)  to  be 


In  resonance  A/T 

Minima  of  Input  resistance  occur  for  kh  given  by 
«  (N  +  )  TT  ,  N  =  1  ,  2,  3,  ... 

min  1 013  ^ 


Cin(2n7T) 
kb  Si(2n7r) 


(3.2-39) 


(3.2-40) 
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37 


The  Input  resistance  at  these  minima  element  lengths  Is  found  from 
Eq.  (3.2-35)  to  be 


( R  ) 

in'^mlnlma 


where 


^  [  y  -  1  +  in[(2N  +  l)7r] 


y =  Euler's  constant  =  0.57721 
The  accuracy  of  Eq.  (3.2-40)  increases  as  N  — ► oo . 

Summary  of  Results 


(3.2-41) 


The  Input  Impedance  and  directive  gain  properties  of  quarterwave  and 
electrlcally-short  monopole  elements  with  groundplanes  of  zero,  large,  and 
infinite  extent  are  compared  in  Table  4  of  Section  3.9.  The  peak 
directivity  is  approximately  3  dB  less  with  groundplanes  of  zero  extent 
than  with  groundplanes  of  large  extent.  However,  the  directive  gain  on  the 
horizon  is  approximately  3  dB  more  with  groundplanes  of  zero  extent  than 
with  groundplanes  of  large  but  finite  extent.  The  radiation  resistance 
with  groundplanes  of  zero  extent  is  approximately  one-half  that  with 
groundplanes  of  large  extent. 

Unlike  dipole  elements  in  free  space  whose  first  resonance  occurs  for 
dipole  half-lengths  approximately  equal  to  a  quarterwave,  monopole  elements 
with  groundplanes  of  zero  extent  have  a  first  resonance  for  an  element 
length  approximately  equal  to  a  half-wave. 


3.3  Integral  Equation,  0  ^  ka  ^  2.75 

In  Section  3.2  the  fields  generated  by  the  monopole  element  were 
determined.  These  fields  Impinge  on  the  groundplane  and  Induce  a 
groundplane  current.  For  sufficiently  small  groundplanes ,  the  fields 
generated  by  the  groundplane  current  may  be  determined  by  Bardeen's 
Integral  equation  method^^^. 

In  this  method  a  cyllndrlcally  symmetrical  electromagnetic  wave 
(generated  by  the  element)  Is  Incident  on  the  groundplane  disc.  The 
fields  generated  by  the  Induced  groundplane  current  are  required.  Th 
solution  depends  upon  solving  an  Integral  equation  of  the  first  kind. 
For  groundplane  radii  of  arbitrary  radius,  the  Integral  equation  Is 
not  readily  solvable  because  It  contains  two  Integrals.  However,  for 
sufficiently  small  groundplane  radii,  Bardeen  neglects  one  of  the 
Integrals  so  that  the  Integral  equation  may  be  solved  explicitly. 

Although  Bardeen  gives  a  general  formulation  of  the  solution  for 
the  resulting  single  Integrand  Integral  equation, his  only  explicit 
results  are  for  the  case  when  the  Incident  wave  Is  generated  by  an 
Infinitely  thin  dipole  element  whose  base  Is  at  a  non-zero  height 
above  the  center  of  the  groundplane.  Bardeen  restricted  his  solution 
to  elements  at  a  non-zero  height  above  the  groundplane  In  order  to 
avoid  having  a  source  point  (the  base  of  the  element)  at  a  near-fleld 
point  of  Interest  which,  for  the  Integral  equation  method.  Includes 
the  center  of  the  groundplane.  The  total  field  Is  then  given  by  the 
vector  sum  of  the  Incident  field  and  the  Induced  field  [see  Eq. 


The  present  case  of  Interest  is  a  tnonopole  element  whose  hase  Is 
In  the  plane  of  the  groundplane,  i.e.,  at  a  vertical  height  v=0  above 
the  center  of  the  groundplane  (see  figure  6).  For  this  case,  the 
total  magnetic  field  may  be  determined  by  first  evaluating  the  field 
with  the  element  at  an  arbitrary  height  v  >  0  and  then  by  evaluating 
the  resulting  expression  In  the  11m  v-*0.  By  such  a  procedure,  an 
Indeterminate  expression  Is  avoided  for  the  field  generated  by  tbe 
Induced  groundplane  current. 


Accordingly,  the  total  magnetic  field  Intensity  f  see  Fq . 

(3. 1-1. 2)]  In  the  limit  v— *0  Is  given  by 

Mm  11m  +  11m  (3.3-1) 

v-*0  V— *0  ^  v-*0 

(e)  (g) 

where  H|j|  and  are  the  magnetic  field  Intensities  generated  by 

the  element  and  groundplane  currents,  respectively.  Tn  the  following 
evaluation  of  the  two  terms  of  Eq.  (3.3-1),  tbe  fields  are  assumed  to 
have  an  e^^*"  time  dependence  [see  Eq.  2.1-2]  unlike  the  time 

dependence  assumed  by  Bardeen. 

The  first  term  of  Eq.  (3.3-1)  is  given  exactly,  for  an  Infinitely 
thin  element,  by  [see  Eq .  3.2-4) 


11m 

v-0 


j  1(0)  r 

47rp8ln(kh)  [ 
i  z  sln(kh) 


exp(-1kr  )  -  cos(kh)  exp(-1kr) 
o 

exp(-ikr)l 


(3.3-2) 


Tn  the  far-fleld,  Eq.  (3.3-2)  reduces  to  [see  Eq.  (3.2-10)] 

lira  =  A  [j  exp(.ikh  cosO)  +  1  cos(kh)  +  cosO  sln(kh)]  , 

V— *0  ^  far-fleld  (3.3-3) 

where 

A  =  - 

o  4  TT  r  sln0  sln(kh) 

The  second  terra  of  Eq.  (3.3-1)  may  be  evaluated  by  utilizing 

(o) 

Bardeen's  Eq.  (31)  for  which  he  obtained  as  a  solution  to  his 

single  Integrand  Integral  equation.  For  an  Infinitely  thin  element 

(e) 

and  for  sufficiently  small  groundplanes  (ka^l),  is  given  by 

Bardeen's  Eq.  (31)  as 


=  sgn(z)  - ^ —  f  F(s)  K(s)  ds,  ka  <  1 

D  0  ( 


where 


(3.3-4) 


-jkA  -jkA 

K(s)  =  e  slnh(kB  )  slnh(ks)  -  e  slnh(kB^)  slnhfka) 


slnh(ka' 


sgn(z)  = 


+1,  z  >  0 
-1,  z  <  0 


e  *  permittivity  of  free  space  (farad/m) 
o 

2  2  1/^  2 

A^  -  jB^  =  fp  +  f-ja  +  z  sgn(z)]  }  =  [r  -  sgn(z)  21ra  cosO  - 

^2  ”  -^®2  °  ^  ^  ^  sgn(z)]^}^^^  =  [r^  -  sgn(z)  21rs  cosO  - 

s  =  dummy  variable  with  the  dimension  of  length 


F(s)  =  excitation  function  related  to  the  radial  electric 

(e)  I 

near-fleld  Intensity  E'  I  „  which  Is  incident  on  the 

r  z^O 

groundplane. 


.  " 


vv 


The  second  term  of  Fq.  (3.3-1)  Is  therefore  given  by 


11m  =  sgn(z) 


=  sgn(z) 


11m  /F(s)K(s)  ds,  ka  <  ] 

•'0 


a 

■/[ 


Urn  F(s)l  K(s)  ds,  ka  1 
V  — *0  J 

(3.3-S) 


The  excitation  function  F(s)  Is  given  by  Bardeen's  Fq.  (3S)  as 

^  f  ^  slnfs  (p^  -  k^)^^^]  .  . 

F(s)  =  -  ^  /  — 7T~TT7T72 - 

f)  (p  -  M 

where 

p  =  dummy  variable  with  the  dimensions  of  (length)  ^ 


h(p)  *  function  of  the  radial  electric  near-fleld  Intensltv 
which  Is  Incident  on  the  groundplane. 

P  Z  “"U 

The  function  h(p)  Is  given  by  Bardeen's  Fq .  (33)  as 


0 


where  J^(x)  =  Bessel  function  of  the  first  kind. 


The  electric  field  Intensity  Incident  on  the  groundplane 

Is  given  by 


»  z-o'  V 


(3.3-B) 


■..  »  X  •  ■  •  •  *  "ji  *■ 


WU>.VJ 


•-  •’> 


where 


=  magnetic  vector  potential  generated  by  the  roonopol e 
source  points 

=  permeability  of  free  space  (henry/m) 

=  unit  vector  In  the  cylindrical  radial  direction 


The  magnetic  vector  potential  A 


(e) 


arbitrary  height  v  Is  given  by 


CI6) 


for  the  element  base  at  an 


where 


.(e) 


u 


■/ 


1(0  (1/PQ)  exp(-jkPQ)  d^ 


/\ 

=  unit  vector  along  the  z  axis 
PQ  =  distance  from  source  point  Q  to  field  point  P 

«  fP^  +  (2  -  ^  ,  Q(x',  y',  O«0f0,  n,  O- 

1(0  =  monopole  current  distribution  *  fl (0) /sln(kh) ]  slnfkfh-^l]. 


In  order  to  facilitate  the  evaluation  of  h(p)  In  Kq .  (3.3-7),  It 

Is  convenient  to  express  the  factor  (1/PO)  exp(-lkPO)  In  terms  of 

^  1711 

Bessel  functions.  Using  Soramerfeld  s  forraula'^  '  and  the  dummy 
variable  p  Introduced  In  Kq.  (3.3-6), 


g-jkPQ 

r  r 

/  ^  ~  '  2  -  ^  -  vl 

4^- 

p  dp 

pq- 

Jq  \  / 

-  (71/2)  <  arg  <  (7r/2), 

-  7r<  arg 

k  <  n 

(1.3-ini 


I 


I.  / 


Substituting  Eq.  (3.3-10)  Into  Eq.  (3.3-9), 

u  00 


1(0)  . 


ATT  sln(kh)  z 


J_(pp)  exp 


(P  -  k^) 


-  (7r/2)  <  arg  Vp  -  k^  ^  (7r/2),  -7r<  arg  k  <  0 


(3.3-11) 


Since  we  are  Interested  In  evaluating  k\  consider  the  case 

p  I  Z“0 

z  <  ^  +  V.  Substituting  Eq.  (3.3-11)  for  such  a  case  Into  Fq.  (3.3-R) 
and  noting  that 


Xe)  -  .  .  S  ’o 

A  u  =  A  and  - - - 

z  z 


A 

,(e)  ^  1  ^ 

>  2=0  -J'^^o^o 


-  pJ^(pp  ), 


where 


-j47ra)fjsin(kh) 


f 

B(P)  =/  sin[k(h  -  ^))  exp 

*'0 


_ 

p^J^(pp)  exp^z-v)  v/p^  - 


B(p)  do 


(3.3-12) 


+  -  k^  sin 


k  exp  [  -hvp  -  k  ]  +  Vp  -  k  sln(kh)  -  k  cos^kh) 

2 

P  (3.3-13) 


Substituting  Eq.  (3.3-13)  Into  Eq.  (3.3-12)  and  letttng  z*0, 

/  -y  r/f\\  /»  T  /  \  '4^  ~  ^  r  -  k 

„(e^  ^  1(0)  /dp  J.(pp)  e  ^  e 

P  1^=0  sln(kh)  I  L 

•^0 


+  -  k^  sln(kh)  -  k  cos(kh)  1 


(3.3-141 


Comparison  of  Eq.  (3.3-14)  with  Eq.  (3.3-7)  yields 

-  (h  +  v)  >4^  ^ 


h(p) 


1(0) 


k  e 


j47rtJ€  sln(kh) 
o 


n.  2 

-  V  Vp  -  k  .  ,  -  V  Vp  -  k 


+  sln(kh)  e  '  ''  "  -  k  cos(kh)  e 


J2  72 1 

vp  -  k 

(3.3-15) 


Substituting  Eq.  (3.3-15)  Into  Eq.  (3.3-6) 
F(s)  = 

where 


-  2  C  k  +  5ln(kh)  I2  -  co8(kh)  T^] 

^  ^  (3.3-16) 


C  s 


1(0) 


1  UTTCJe  sln(kh) 


Jo 


sln(  s  -  k^)  dp 


2  .  2\  1/2 


(p^- 


/ 


00 


p  sln(s 


-2.  .  „ 


4^  -  k^)  g-  V  Vp  -  k 


/ 


00 


p  sln(s  4^  “  1^^)  e 
2  , 2\  1/2 


sP-  72 

-  V  Vp  -  k 


dp 


(p^- 


In  Integrals  I2,  I^,  we  Introduce  a  change  In  variable  p  to  y  given  by 
y  -  (p^  -  k^)  -  (7r/2)  <  arg  y  <  (7r/2)  (3.3-17) 


where  the  condition  on  y  follows  from  the  conditions  of  Eq .  (3.3-10). 


46 


From  Integral  2.663.1  of  Ref.  [22],  it  follows  that  Is  given  by 


00 

At. 


8ln(sy)  dy  =  -- sin(sy) 

(h  +  v)  +  s 


-  s  cos(sy)] 


-k(h  +  v) 

- 9 - 5“  [~j(h  +  v)  slnh(ks)  -  s  cosh(ks)] 

(v  +  h)  +  s 

(3.3-lR) 


From  Integral  2.667.5  of  Ref.  [22],  It  follows  that  is  given  by 

_  -vy  r  /  2  2  \ 

y  sln(sy)  e  dy  =  ^  (-vy  +  ~ —  j  slnfsy) 

v+s  [\  v+s'  J 


sy  +  -j - Y'  )  <^os(sy) 

8  +  V  / 


"jlcv  2  2 

-  kv  Slnh(k8)  -  Jks  coshCks)  -  -  .«■  >  syfTcs)  -  2  vs  coshCks-> 

3  +  V  S^+  V^ 


(3.3-1*)) 


From  Integral  2.663.1  of  Ref.  [22]  It  follows  that  I^  Is  given  by 

/°° 

8ln(sy)  e'^^y  dy  =  e"^’^  _  v.A^n(sy)  -  s  cos(sy) 

s  +  V 


(3.3-20) 


-jkv  ~  j  V  sinh(ks)  -  s  cosh(ks) 

2.  2 
S  +  V 


Consider  now  the  limits  of  I^,  I^,  and  I^j  as  v— >0. 

-Ikh 

lira  I.  =  - 5 - ^ -  [-J  h  slnh(ks)  -  s  cosh(ks)]  . 

v-0  h  +  s  ( 
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11m  I_  =■  -7rk(5(8)  slnh(ks)  + 
v-0 


+  jTT  <5*^(8)  slnh(ks)  - 


11m  I-  =  j7T5(s)  slnh(ks)  + 

v-0 


j  k  co8h(ks) 


J  slnh(ks) 


cosh(ks) 


-7r5'(8)  C08h(k8) 

r'1.3-22’) 


(3.3-23) 


where 


r  1  V 

<5(x)  •  Dirac  delta  function  =  11m  j'-'^  ~2 - 2~ 

V— *oL  V  +  X 


d(x-x  )  dx 
o 


^ '(x) = 


1. 

a  <  X 

o 

<  b 

1/2, 

a  =  X 

o 

<  b  or  a  < 

.  0, 

X  <  a 
o 

or  X  >  b 
o 

lim 

a 

1  V 

v-0 

^  v^H“  X 

In  Eq8.  (3.3-22)  and  (3.3-23),  tho8e  terms  containing  the  product 
d(s)  8lnh(ks)  may  be  8et  equal  to  zero  eince  from  Eq.  (B.  2-12)  of  Ref. 
[23],  X  ^(x)  »  0. 


The  11m  F(8)  Is  found  by  substituting  Eqs.  (3.3-21)-(3.3-23)  Into  '^q. 
v—O 


(3.3-16). 

Accordingly, 


11m  F(8) 
v-0 


2Ck 

^-jkh 

f-.1  h  slnh(k 

n 

2 

h  +  8 

2C 

sin(kh) 

-j  k  cosh  (ks) 

n 

s 

2C 

k  cos(kh) 

cosh(ks) 

2 

s 


(3.3-24) 


The  magnetic  field  Intensity  generated  by  the  groundplane  current,  for 
the  case  v-*0.  Is  found  by  substituting  Eq.  (3.3-24)  Into  Eq.  (3.3-S).  The 
total  magnetic  field  Is  explicitly  determined  by  numerical  evaluation  of  Eqs. 
(3.3-1)  -  (3.3-5). 

In  the  far-fleld,  the  factor  K(8)  and  the  parameters  ,  E, ,  In 

Eqs.  (3.3-4)  and  (3.3-5)  reduce  to  a  simpler  form. 

In  the  far-fleld  when  r— »»,  then  a  <<  r  and  s  <<  r.  Eor  these 
conditions , 

2  2  2  1/2 
-  jBj^~[r  -  sgn(z)  2jra  cos  ©  -  a  cos  ©1  =  r  -  sgn(z)  Ja  cos©, 

a  «  r  (3.3-25) 

^2  ”  ~  S8n(z)  2jra  cos  ©  -  s^cos^©]^^^  =  r  -  8gn(z)  ,1s  cos©, 

s  «  r  (3.3-26) 

Equating  real  and  Imaginary  parts  In  Eqs.  (3.3-25)  and  (3.3-26) 

A^  -  Aj  »  r;  a  «  r,  s  «  r  (3.3-27) 

=  sgn(z)  a  cos  ©,  a  «p  (3.3-2R) 

^2  “  sgn(z)  s  cos  ©,  s  «p  (3.3-?©) 

Substituting  Eqs.  (3.3-27),  (3.3-28),  (3.3-2©)  Into  Fqs.  (3.3-4)  and  (3.3-5), 
the  far-fleld  magnetic  field  Intensity  due  to  the  groundplane  cxirrent,  for 
the  case  v— *0,  la  given  by 


exp(-jkr) 
kr  sin© 


11m  F(s) 
v-0 


g(s)  ds;  far-fleld,  ka  ^  1 
(3.3-30) 


llm 

v-O  ♦ 


where 


g(s) 


slnh(ka  cosQ)  slnh(ks)  -  slnh(ks  cosQ)  slnhfka') 
sinh(ka) 


and  lim  F(s)  Is  given  by  Kq .  (3.3-24) 
v-0 


Eq .  (3.3-30)  can  be  reduced  further  by  utilizing  the  properties  of  the 
Dirac  delta  function  given  by  Eq .  (R.2-7)  of  kef.  (23]  and  the  relation 


J'  d'(x)  f(x)  dx  - - ^  f'(0'''). 


It  follows  that  Eq .  (3.3-30)  reduces  to 
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(3.3-31) 


where  A^  and  g(s)  are  defined  In  Eq.  (3.3-3)  and  Eq.  (3.3-30), 
respectively. 


The  total  far-fleld  magnetic  field  intensity  generated  by  the 
element  and  groundplane  currents,  for  the  case  v— »0,  Is  found  by 
substituting  Eq.  (3.3-3)  and  Eq.  (3. 3. -31)  Into  Eq.  (3.3-1).  The 

resulting  expression  for  lira  -jg  if,  a  form  suitable  for 


computer  evaluation. 


The  radiated  power  density  s(9),  directive  gain  d(9)  ,  and 
radiation  resistance  R  are  found  from  Eqs.  (2.2-4)  and  (2.2-7).  These 
quantities  are  computed  In  program  "BARDEEN"  in  Appendix  B1 .  Computer 
printouts  of  the  relative  power  radiation  pattern  radiation 

resistance  R,  and  the  directive  gain  d(9)  are  given  In  tables  (Al-l  - 
Al-13),  (Al-14  -  Al-25),  (Al-26  -  Al-37),  (Al-38  -  Al-49)  for 
normalized  element  lengths  h/A  =  0.25,  0.1,  0.040,  0.01,  respectively, 
and  normalized  groundplane  radii  =  0,  0.25,  0.50,  ...  3.0  wavenumbers. 

The  radiation  patterns  have  no  appreciable  change  In  shape  for 
groundplane  radii  0  £  ka<1.75  and  resemble  that  of  a  dipole  in  free 
space  with  peak  gains  approximately  In  the  direction  of  the  horizon 
and  with  directivities  less  than  that  of  a  dipole  whose  total  element 
length  is  twice  that  of  the  monopole  element  length.  For  a 
quarterwave  monopole  element,  the  directive  gain  on  the  horizon 
deceases  from  1.88  dBl  for  ka  ••  0  to  1.23  dBl  for  ka  =  1.75, 

In  table  5  of  Section  3.9,  the  radiation  resistances  obtained  by 
the  integral  equation  method  are  compared  with  those  obtained  by  the 
method  of  moments  for  a  quarterwave  element  and  groundplane  radii 
0  _<  ka  £  3.0.  The  values,  obtained  by  the  Integral  equation  method, 
differ  from  those  obtained  by  the  method  of  moments  by  less  than 
1%  for  0  £  ka  ^  1.75  and  by  less  than  lOX  for  0  £  ka  £  2.75.  These 
results  suggest  that  the  Integral  equation  method  is  accurate  for 
groundplane  radii  0  £  ka  £  1.75  and  is  useful  for  0  ^  ka  £  2.75.  It 
Is  suspected  that  the  best  available  results  are  obtained  by  the 
integral  equation  method  for  0  _<  ka  £  1.0. 
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3.4  Method  of  Moments,  0  <  ka  <14 


In  Section  3.3  it  was  shown  that  the  integral 

equation  method,  for  determining  the  fields  generated  by  the 
groundplane  current  when  the  current  distribution  on  the  monopole 
element  is  initially  specified.  Is  accurate  only  for  groundplane  radii 
less  than  approximately  1.75  wavenumbers.  In  Section 
3.5  it  will  be  shown  that  the  algorithms  utilized  in  the  oblate 
spheroidal  wave  function  metho^t  are  accurate  only  for  groundplane 
radii  no  smaller  than  3.0  wavenumbers.  The  question  arises:  What 
method  is  accurate  over  a  range  of  groundplane  radii  which  includes 
the  region  1.75  <  ka  <  3.0? 

We  have  found  that  Richmond's  method  of  moments^^^  is  the  only 
method  In  the  present  literature  to  be  accurate  over  a  range  of 
groundplane  radii  which  Includes  the  region  1.75  <  ka  <  3.0.  Although 
this  method  Is  primarily  Intended  for  use  when  the  current 
distribution  on  the  monopole  element  Is  Initially  unknown,  this  method 
Is  also  applicable  when  the  element  Is  specified  to  have  a  sinusoidal 
current  distribution. 

This  method  Is  discussed  In  Section  4.2.  In  this  method  the 
element  is  subdivided  Into  N  equal  segments  and  the  groundplane  Is 
subdivided  Into  M  concentric  annular  zones  of  equal  width.  The 
unknown  current  distributions  on  the  element  and  groundplane  are 
expanded  as  a  series  of  N  +  M  overlapping  sinusoidal  dipole  modes 
( slnusoldal-Galerkln  method)  each  with  an  unknown  current 
distribution.  The  N  +  M  currents  are  determined  by  inversion  of  a 
(N  +  M)  X  (N  +  M)  matrix.  The  numbers  of  subdivisions,  N  and  M,  are 


limited  by  the  cost  of  computation  time  and  by  the  precision  of  the 
computer.  The  accuracy  of  the  solution  can  be  decreased  appreciably 
If  either  N  or  M  Is  too  small  or  too  large.  The  method  of  moments 
converges  to  a  solution  when  an  Increase  or  decrease  of  unity  In  the 
value  of  N  or  M  does  not  appreciably  alter  the  solution  for  Input 
Impedance. 

A  particularly  useful  property  of  the  sinusoidal  -  Galerkln  method 
Is  the  sinusoidal  current  distribution  which  Is  Imposed  on  the  element 
by  setting  N=l. 

The  Input  Impedance  and  radiation  pattern  of  thin  quarterwave 
elements  (b/X  =  10  h/ X  =■  0.25),  for  groundplane  radii  0  <  ka  14 
wave  numbers,  were  determined  by  Richmond's  method  of  moments  by 
utilizing  MITRE  programs  "RICHMDl”  and  ”RICHMD2"  which  are  discussed 
In  Section  4.2  and  listed  In  Appendix  B5.  With  N=l,  convergent 
solutions  were  obtained  for  values  of  M  given  by: 

M  =  3,  ka  »  0.25 

M  -  7,  ka  =  0.50 

M  =  16,  ka  -  0.75,  1.0  .  8.5 

M  =  3ka,  ka=*8.75,  9.0,  .  14.0. 

Computer  printouts  of  the  Input  admittance,  groundplane  current 
distribution,  radiation  resistance,  directive  gain,  and  radiation 
patterns  are  given  In  Appendix  A5. 

The  Input  resistance  and  reactance  —  as  determined  by  the  method 
of  moments  or  methods  which  give  similar  results  —  are  plotted  for 
thin  quarterwave  elements  on  groundplanes  of  radii  0  £  ka  £  14  In 


figures  9  and  10  of  Section  3.9.  The  directive  gain  on  the  horizon, 
peak  directivity,  and  elevation  angle  of  the  peak  directivity  given  by 
the  computer  printouts  of  Appendix  A5  are  plotted  In  figures  11  -  13 
of  Section  3.9. 

The  Input  Impedance  of  thin  elements  with  a  sinusoidal  current 
distribution  were  also  determined  for  element  lengths  h/A  =  0.1  and 
0.025  and  groundplane  radii  0  <  ka  <  8.0  wavenumbers  by  utilizing 
program  RICHMDl  with  N=l.  The  radiation  resistance  of  these  elements 
is  compared  In  figure  14  of  Section  3.9  with  those  of  a  quarterwave 
element.  In  figure  14  the  radiation  resistance  Is  normalized  to  the 
value  of  radiation  resistance  of  each  element  for  ka  =  0. 

Values  of  the  radiation  resistance  for  quarterwave  elements  are 
determined  In  Appendix  A5  by  matrix  inversion  (program  RICHMDl)  and 
also  by  the  far-fleld  radiation  pattern  (program  RICHMD2).  The  values 
determined  by  both  methods  differ  by  less  than  1%  for  small 
groundplane  radii  and  differ  by  less  than  3?  for  the  larger 
groundplane  radii. 

The  values  of  radiation  resistance  obtained  from  program  RICHMD2 
are  compared  in  table  5  of  Section  3.9  with  the  values  obtained  by  the 
Integral  equation  method  and  the  oblate  spheroidal  wave  function 
method.  Richmond's  method  of  moments  gives  useful  results  over  the 
entire  range  0  <  ka  £  14  and  gives  good  agreement  with  the  Integral 
equation  method  for  0  <  ka  ^  1.75  and  with  the  oblate  spheroidal  wave 
function  method  for  3.0  ^  ka  ^  6.5.  Whereas  the  method  of  moments 
gives  useful  results  In  the  regions  1.75  <  ka  <  3.0  and  6.5  <  ka  £  14, 
the  other  methods  fall  in  these  regions.  For  ka  >  14,  Richmond's 
method  of  moments  Is  not  as  useful  because  of  Increased  computation 
time  and  decreased  accuracy.  It  is  suspected  that  for  thin  monopole 


elements  the  best  available  results  are  obtained  by  the  method  of 
moments  for  1.25  £  ka  ^  2.75  and  6.75  £  ka  £  14.  For  relatively  thick 
monopole  elements,  Richmond's  method  of  moments  gives  the  best 
available  results  for  groundplane  radii  0  <  ka  <  14  as  discussed  In 
Section  4.2. 

3.5  Oblate  Spheroidal  Wave  Functions,  3.0  £  ka  £  6.5 

Oblate  spheroidal  coordinates  (see  table  1  In  Section  2.1)  are 
particularly  convenient  for  handling  the  boundary  conditions  of  the 
magnetic  field  Intensity  on  the  groundplane.  The  requirement  for 
constant  tangential  magnetic  field  Intensity  across  the  upper  and 
lower  hemispheres  at  the  groundplane  interface  may  be  specified  at  all 
points  of  the  groundplane  disc  by  a  boundary  cotidltlon  at  only  the 
oblate  "radial"  coordinate  ^  =  0. 

For  groundplane  radii  of  the  same  order  of  magnitude  as  the 
excitation  wavelength,  Leltner  and  Spence^^^  utilized  oblate 
spheroidal  wave  functions  to  determine  the  groundplane  current  Induced 
by  a  thin  quarterwave  element  with  a  sinusoidal  current  distribution. 
Leltner  and  Spence  give  numerical  values  of  the  groundplane  current 
distribution,  radiation  resistance,  and  far-Fleld  power  density  (at 
constant  element  base  current)  for  groundplane  radii  ka  =*  3,  4,  5 ,  ^4^ 
wavenumbers.  The  complex  current  distributions  on  both  the  top  and 
bottom  faces  of  the  groundplane  are  reported. 

In  this  section,  we  report  the  results  of  a  computer  program 
"MONOPL,"  based  on  the  theory  of  Leltner  and  Spence,  which  calculates 
the  directivity  pattern  and  radiation  resistance  for  groundplane  radii 
Including  the  cases  ka  =  3,  4,  5,  wavenumbers.  Our  results  are 
consistent  with  (but  not  Identical  to)  the  results  reported  by  Leltner 
and  Spence.  The  form  of  the  solution,  corrections,  regions  of 
calculation  validity,  accuracy,  and  numerical  results  of  program 
MONOPL  are  discussed  In  the  remainder  of  this  section. 


Form  of  Solution 


Substitution  of  Kqs.  (3.5-1)  and  (3.5-2)  Into  Eq.  (2.2-4)  yields  the 
directivity  d'(?7)  given  by 


«'(«  -  2  I  i„  r 


(3.5-3) 


Corrections 


Problems  were  encountered  in  generating  the  eigenvalues  of  the 

oblate  spheroidal  wave  functions.  These  problems  were  narrowed  down 
to  the  continued  fractions  that  generate  the  eigenvalues.  These 
continued  fractions  were  derived  from  the  recurrence  relations  which 
In  turn  come  from  the  differential  equation.  A  discrepancy  was 
discovered  between  our  derived  continued  fractions  and  the  published 
continued  fractions  [3]  .  The  corrected  continued  fractions  appear 
below. 
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The  eigenvalues  can  also  be  expressed  In  terms  of  a  series  expansion, 


which  has  the  form 


°  Z  ^2k  ^ 

k=l 


£m  2k 


(3.5-6) 


The  first  two  terms  In  this  expansion  were  checked  against  the  continued 

fractions  and  were  found  to  agree.  This  Is  Important,  since  the  continued 

fraction  method  by  which  values  of  the  eigenvalues  are  obtained  depends  upon 

the  accuracy  of  the  roots  V,  In  equation  (3.5-6). 
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Regions  of  Calculation  Validity 

A  lower  bound  on  the  value  of  f  =  2TTa/X  for  which  the  calculations  are 
valid  depends  upon  the  following  equations  [see  Kqs.  (34)  and  (4R)  of 
Ref.  31; 
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where  the  prime  indicates  summation  over  alternate  n,  starting  with  n=0  if 
(i  -  1)  Is  even,  with  n=l  If  (  i  -  1)  Is  odd. 

Eq.  (3.5-7)  Is  the  expression  for  the  first  order  radial  spheroidal 
wave  function  of  the  first  kind.  Eq.  (3.5-8)  is  used  in  another 
expression  to  obtain  the  first  order  radial  spheroidal  wave  function 
of  the  second  kind.  As  €  becomes  small,  ^  becomes  large,  so  more 
terms  are  needed  in  Eqs .  (3.5-7)  and  (3.5-8).  In  theory,  Eq.  (3.5-7) 
and  (3.5-8)  converge  for  all  real  values  of^*7r/2€.  Computationally, 
however,  because  of  the  finite  accuracy  of  the  computer,  Eqs.  (3.5-7) 
and  (3.5-8)  will  not  converge  for  all  real  values  of 
^  =  TT/2€.  To  be  on  the  safe  side,  we  can  restrict  ^  so  that  ^  Is  less 
than  unity.  With  this  assumption,  we  get  the  computational  constraint 
^  =  7Tl2€  <  1  which  implies  a  lower  bound  on  e  given  by  6  >  n’/2  w  1 . 57 . 

An  even  tighter  lower  bound  on  €  Is  obtained  by  observing  what 
happens  In  the  algorithm  that  is  used  to  obtain  the  eigenvalues.  This 
algorithm  Is  not  well  behaved  for  values  of  €  less  than  €  =  2.5. 

Therefore  e  Is  lower  bounded  by  €  2.5.  For  accurate  values  of 

radiation  resistance,  i  is  lower  bounded  by  €  3.0  (see  table  B  of 

Section  3.9). 


An  upper  bound  on  6  Is  obtained  by  observing  what  happens  in  the 
continued  fraction  algorithm  that  is  used  to  obtain  the  eigenvalues. 
For  values  of  (.  greater  than  6.5,  the  series  expansion  for  the 
eigenvalues,  given  by  Eq.  (3.5-6),  does  not  give  an  accurate  enough 
answer  using  only  the  first  four  tenns.  The  resulting  values  for  the 
eigenvalue  are  far  enough  from  the  correct  eigenvalues  to  cause  the 
continued  fraction  algorithm  of  Eqs.  (3.5-4)  and  (3.5-5)  to  converge 
to  a  root  which  Is  not  the  eigenvalue.  As  a  consequence,  erroneous 
values  of  radiation  resistance  can  be  obtained  for  € >  6.5  (see  table  5 
of  Section  3.9). 

Consequently,  the  range  of  £  for  which  Leltner  and  Spence's  method 
of  oblate  spheroidal  wave  functions  Is  useful  is  2.5  ^  ^  <  6.5.  It  Is 
suspected  that  the  best  available  results  are  obtained  by  the  method 
of  oblate  spheroidal  wave  functions  for  3.0  <  €  ^  6.5. 

Accuracy 

There  are  many  equations  involved  in  the  calculation  of 
directivity  and  radiation  resistance.  Some  of  these  equations  involve 
series  expansions.  When  we  varied  the  number  of  terms  In  these 
series,  the  radiation  resistance  was  found  to  vary  only  In  the  fifth 
or  sixth  significant  figure. 

Another  problem  that  was  mentioned  previously  was  the  accuracy  by 

which  Eq.  (3.5-6)  computes  the  eigenvalues.  Because  of  computational 

reasons,  only  the  first  four  terms  In  Eq.  (3.5-6)  were  used.  Because 

of  this,  a  raw  eigenvalue  is  computed  using  Eq.  (3.5-6),  which  Is  then 

used  as  an  Initial  guess  In  the  continued  fractions.  The  continued 

fractions  have  many  roots  In  V  •  The  number  of  roots  Is  dependent 
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upon  the  number  of  fraction  terras  used  In  the  continued  fraction. 

Only  one  of  these  roots  however  can  be  the  eigenvalue.  If  the  raw 
eigenvalue  Is  far  enough  from  the  correct  eigenvalue  then  the 
continued  fraction  will  converge  to  a  root  which  Is  not  an  eigenvalue. 
This  will  produce  wrong  results. 


The  eigenvalues  that  were  computed  were  checked  against  published 
values^“^^-  The  computed  values  were  found  to  be  within  the 
percentage  error  of  the  published  values. 

The  computed  values  of  directivity  on  the  horizon  and  radiation 
resistance,  for  different  values  of  groundplane  extent,  were  found  to 
agree  with  the  relationship  given  by  Eq.  (2.3-3)  to  at  least  five 
places  after  the  decimal  point  for  the  free  space  wave  Impedance 
^  =  376.73037  ohms. 

Numerical  Results 

The  far-fleld  power  density  s'(?7),  radiation  resistance  R,  and 
directive  gain  d'(T])  of  quarterwave  elements  given  by  Eqs.  (3.5-1), 
(3.5-2)  and  (3.5-3),  respectively,  were  numerically  evaluated  by 
program  MONOPL  written  In  FORTRAN  77  language  for  use  on  a  DEC 
PDP-11/70  computer.  The  program  listing  Is  given  In  Appendix  B3. 
Numerical  values  were  obtained  for  the  cases  2.5  £  ka 

The  computed  eigenvalues  are  given  In  table  2.  Computer  printouts 
of  the  directivity  patterns  are  given  In  Appendix  A3  for  ka  =  3,  4, 

5,  y42 .  The  patterns  are  plotted  In  figure  8  of  Section  3.9  as  polar 
graphs  on  the  same  linear  scale.  In  these  plots  the  total  radiated 
power  Is  held  constant. 

The  directional  gain  on  the  horizon,  radiation  resistance  and  peak 
directivity  are  summarized  In  table  3. 
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3.6  Scalar  Diffraction  Theory,  Geometric  Theory  of  Diffraction, 


6.5  <  ka  <  00 


For  an  element  on  a  groundplane  of  sufficiently  large  radius, 
Tang^^^  utilized  a  scalar  theory  of  diffraction  to  calculate  the 
far-field  elevation  pattern.  For  elevation  angles  near  the  horizon 
(9a7r/2)  the  far-fleld  magnetic  field  Intensity  is  determined  by 
linear  extrapolation  to  the  result  for  the  element  itself  given  by 
Eq.  (3.2-10)  with  0*  7r/2  rad.  Tang's  method  for  the  radiation  pattern 


Is  more  accurate  than  that  obtained  by  the  variational  method  of 


Section  3.7  because  it  Includes  an  additional  term  In  the  expansion 
for  the  total  magnetic  field  intensity.  Since  the  variational  method 
is  useful  for  groundplane  radii  as  small  as  ka=30  wavenumbers,  Tang's 
method  should  be  useful  for  even  smaller  groundplane  radii  provided 


ka  »  1. 


The  geometric  theory  of  diffraction  (GTD)  Is  another  method  which 
is  applicable  for  sufficiently  large  groundplane  radii.  In  GTD,  the 
fields  are  determined  by  ray  optics  (an  Incident  ray  plus  a  reflected 
ray)  and  diffraction  by  the  edge.  However,  the  effect  of  edge 
diffraction  Is  only  approximated  because  In  GTD  the  edge  Is  treated 
point  by  point  as  though  it  were  a  straight  knife  edge  of  Infinite 
extent.  For  this  reason,  GTD  may  be  applied  to  an  element  at  the 
center  of  a  circular  groundplane  only  when  the  groundplane  Is  of 
sufficiently  large  radius.  The  method  of  GTD  is  reviewed  by 


Balanls''  ■'  who  also  gives  a  computer  program  for  calculating  the 


diffraction  coefficient. 


In  Section  4.3,  the  method  of  moments  combined  with  GTD  yields 
results  for  input  Impedance  which  are  useful  for  groundplanes  of  radii 
ka  >  6  wavenumbers  and  are  accurate  for  ka  >  8.  Therefore,  when  the 


element  current  distribution  is  constrained  to  be  sinusoidal,  the 


method  of  GTD  Is  expected  to  give  useful  results  for  groundplane  radii 
ka  >  6.5  wavenumbers  over  which  range  the  method  of  oblate  spheroidal 
wave  functions  does  not  give  useful  results. 


3.7  Variational  Method,  30  <  ka  <  oo 


For  an  element  on  a  groundplane  of  very  large  radius 
(ka  >  30  wavenumbers),  differences  In  the  Input  impedance  and 
radiation  pattern  from  that  for  a  groundplane  of  infinite  extent  may 
be  determined  by  utilizing  a  variational  method  of  Storer^^^ ’ 

With  reference  to  figure  1  of  Section  2.1  and  for  a  sinusoidal 
current  distribution  on  the  element  given  by  Fq .  (3.1-1),  the  Input 
Impedance  difference  Is  given  by  Eq .  (20)  of  [7)  as 


Z 


in 


Z  » 
00 


where 


-j  V exp( j2ka) 
ATTka 


1  -  cos(kh) 
sln(kh) 


j  ^  exp[i(2ka  + 
(ATTka) 

ka  >  30  (3.7-1) 


Z^^  =•  Input  Impedance  for  an  element  on  a  groundplane  of 
radius  ka  wavenumbers  (ohms) 

Z^*  Input  impedance  for  the  same  element  on  a  groundplane 
of  radius  ka  =oo wavenumbers  (ohms) 

^  =  wave  Impedance  In  free  space  =  376.73037  ohms 

1  /2 

Since  (ATTka)  »  1  for  ka  ^  30,  the  Input  resistance  difference 

R.  -  R_^  and  Input  reactance  difference  X,  -  X  are  given 
In  00  ^  In  00 

approximately  by 


R 


in 


'  00 


X 


In 


X 


00 


si n(2ka) 
ATTka 

cos(?ka) 

ATTka 


ka  >  30 

(3.7-2) 


ka  >  30 

(3.7-3) 


37r/A)] 
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The  maxima,  nulls,  and  minima  of  —  R  occur  for  values  of  ka  given 
approximately  by 


ka  = 


(2N  +  ^  )  (7r/2) 
-  N  (7r/2) 

(2N  +  ^  )(7r/2) 


(maxima) 

(nulls)  ,  N  =  0,1,2, 
(minima) 

(3.7-4) 


The  Input  Impedance  given  by  Eq.  (3.7-1)  lo  calculated  In 

computer  program  "MONOSTOR"  whose  listing  Is  given  In  Appendix  B-4. 


For  very  thin  quarterwave  elements,  36.54  ohms  and 

21.26  ohms  (see  Section  3.8).  The  numeric  directivity  on  the 

horizon,  d(7r/2).  Is  related  to  the  radiation  resistance  R  =  R,  by 

In 

Eq.  (2.3-2).  Computer  printouts  of  the  maxima  and  minima  of  R,  -  R  . 

in  OCr 

Xin  “  given  In  Appendix  A4.  For  very  thin  quarterwave 

elements.  Appendix  A5  also  gives  computer  printouts  of 
d(7r/2),  and  D(rr/2)  -  10  log  d(7r/2). 

Differences  In  the  far-fleld  radiation  pattern  from  that  for  a 
groundplane  of  Infinite  extent  are  given  In  [8].  For  the  waveform  of 
Eq .  (2.1-2)  and  an  element  sinusoidal  current  distribution  I(z)  (see 
figure  I  of  Section  2.1),  the  difference  In  the  far-fleld  magnetic 
field  Intensity  Is  given  by  Eq.  (6)  of  (8]  as 


. 


mkV  »  - 
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“  ”(ll^ka=oo° 


sgn  e  r  j  ^(0)  exp(-.ikr)  [l-cos(kh)1 
L  sln(kh)  27rr  ■  TH 


exp[-jka(l-stn9  co8(|l)1  co8(l)  d(ll 


(l-sln0  cos(tl) 


J  .  ka  > 
(3.7-5) 


where 


Hjji  =  far-field  magnetic  field  Intensity  for  an  element  on  a 
groundplane  of  radius  ka  wavenumbers  (amp/m) 

=  far-fleld  magnetic  field  Intensity  for  an  element  on  a 

00 

groundplane  of  Infinite  extent  (amp/m) 


sgn  6 


+1,  0  <  0  <  7r/2 


-1 ,  rr/2  <  0  TT 


The  magnetic  field  Intensity  ^(^Ii^^g^ools  given  by  (see  Section  3.R) 


^(>  ^  ka=  00 


i  1(0)  exp(-.1kr) 
sln(kh)  27rr 

0,  7r/2  <  0  <  TT 


C03(kh  COS0)  -  co8(kh) 
sin  0 


The  far-fleld  electric  field  Intensity  =•  where  H|j|  Is  given  by 
Eq  (3.7-5)  and  V  Is  the  free-space  wave  Impedance. 


One  of  the  distinguishing  features  of  the  far-fleld  radiation 
pattern  for  groundplane  of  large  radii  is  the  occurrence  of  a 
fine-structured  loblng  pattern,  H|j|  --  superimposed  on  the 

pattern  for  a  groundplane  of  Infinite  extent.  The  loblng  pattern, 
given  by  Eq.  (3.7-5),  Is  s)rmmetrlcal  about  the  horizon  (expect  for  a 
phase  shift  ofTTrad)  with  the  most  prominent  lobes  near  the  zenith  and 
nadir  directions.  The  n^^  maximum  of  the  loblng  pattern  in  each 
quardrant  decreases  with  Increasing  values  of  n  where  n«l  corresponds 
to  the  lobe  nearest  the  zenith  or  nadir  direction. 


I 


iJJJl 

I 


w 


The  lobing  structure  Is  characterized  in  [8]  by  the  elevation 
angle  0^  of  the  n^^  maxima,  the  angular  separation  40  between  lobes, 
the  angle  within  which  all  prominent  maxima  occur,  and  the  number 

N  of  prominent  lobes.  The  elevation  angle  0^  in  radians  is  given 
approximately  by 


0.59  (7r/ka),  n=l 

,  ka  >  30 

(n  +  0.75)(7r/ka),  n  =  2,  3 . N 


(3.7-7) 


The  angular  separation  4®  in  radians  is  given  approximately  by 


1.16  (TT/ka),  n=l 


7T /ka,  n  =  2,  3,  . . .  .N 


,  ka  >  30 


(3.7-8) 


Prominent  lobes  are  defined  in  f8]  as  lobes  whose  maxima  are  less  than 
one-fourth  the  amplitude  for  that  of  a  groundplane  of  infinite  extent. 
The  angle  in  radians  is  given  approximately  by 


9  »  1.87  r  *  ~ 

max  L  kh  sln(kb) 


-,2/3 

,  ka  >  30 

-I  ~  (3.7-Q) 


The  number  N  of  prominent  lobes  is  determined  by 


(1.7-10) 


where  6,, 40,  and  0  are  given  by  Kqs.  (3.7-7)  -  (3.7-9),  respectively. 
''1'  ’  max  o  y  -1 


Solving  for  N, 


N  -  °  1.87(7r/ka)  '  -0.25,  ka  >  30 


(3.7-11) 


m 


m 


3.8  Method  of  Images,  ka  ■  oo 


For  the  Idealized  case  of  a  monopole  element  mounted  on  a  groundplane 

of  Infinite  extent  and  of  Infinite  conductivity,  the  monopole  antenna  may 

be  modelled  by  the  method  of  Images  as  a  dipole  with  one-half  the  Input 

f  ?  SI 

Impedance  and  double  the  peak  numeric  directivity  of  the  dipole.  "Hie 

Infinite  groundplane  prevents  monopole  radiation  Into  the  hemisphere  below 
the  groundplane  but  allows  a  radiation  pattern  identical  to  that  of  the 
dipole  In  the  upper  hemisphere. 


In  this  section.  It  Is  assumed  that  the  current  has  a  waveform  given  by 
Eq .  (2.1-2)  and  a  current  distribution  I(z)  on  the  element  and  Its  Image 
(see  figure  7)  given  by 


T(z) 


[l(0)/sln(kh)]  slnfk(h-z)],  0  ^  z  ^  h  (element) 

(3.8-1) 

[l(0)/aln(kh)]  sln[k(h+z)],  -  h  <  z  <  0  (Image) 


The  near-flelds,  far-flelds,  and  Input  Impedance  of  an  electrlcall y-tb1 n 
element  on  a  groundplane  of  Infinite-extent  are  summartzed  In  the  remainder 
of  this  section. 

Near-fields 


The  exact  magnetic  field  Intensity  H  =  and  electric  field 

intensity  E  =  Infinitely  thin  element  are  given  with 

reference  to  Fqs.  (7-15)  and  (7-17)  of  r251,  as 


H 


j  1(0)  exp(-ikR  )  +  exp(-1kR  )  -  2  cos(kh)  exp(-ikr) 
47rp8ln(kh)  ^  . 

^  7.  y  {} 


0,  z  <  0 


(3.8-2) 
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■1  ^?T(0) 
47rpsln(kh) 


(z+h)  exp(-1kR  ) 


where 


= 


F,  = 


f-  (z-h)  exp(-lkR^) 


L.  I  y 

-  2z  cos(kh)  exp(-lkr)  |  ^  ^ 

T  J  ’  ^  ^  " 


0,  z  <  0 

1(0) 

ATT  sln(kh) 


exp(-ikR^)  expC-.ikR^) 


R, 


1  2 

-  2.  cos(kh)  exp(-  ikr) 
r 

0,  z  <  0 


>  0 


R,  =  fp  +  (z-h)  ] 


2,1/2 


R2  =  TP  +  (z+h)  1 


2,1/2 


(3.R-3) 


(3.R-4) 


The  magnetic  field  Intensity  at  the  top  and  bottom  surfaces  of  the 
groundplane  are  given  by  Eq .  (3.8-2)  as 


”(t'z-0^  °  Wrsln(kh)[®’'P^-^’"'^3^  ■  exp(-1kp) 

top  of  groundplane 

^()l^z=0  ~  groundplane 


where 


I!,  -  +  hb'« 


(3.8-5) 

(3.8-6) 


The  Incremental  groundplane  current  4I(P)  contained  within  a  differential 
azimuthal  angle  d(tl  may  be  de  ermlned  from  Fqs.  (3.8-5),  (3.8-6),  and 
Ampere's  circuital  law  applied  to  a  closed  path  along  the  top  and  bottom 
surfaces  of  the  groundplane  along  arc  lengths  within  d(|i.  Since  the 
groundplane  current  I(P)  Is  defined  to  be  positive  In  the  posl 1 1 ve  P dl recti  on 
(see  figure  1  of  Section  2.1),  the  path  is  taken  In  the  clockwise  direction. 


Accordingly, 


41(P)  -  (-u^  ■  p  .fit  -  f-H4,l^.„  +  "^l^.„J  P  H4I 

('3.R-7) 


The  total  groundplane  current  T(p)  Is  given  by 


2  7T 


I(P>  -  f  4I(P)  d*  - 


J  1(0)  fcos(kh)  exp(-lkp)  -  exp(-ikR  )1  (l.R-R) 


sln(kh) 


Eq.  (3.8-ft)  agrees  with  that  given  by  Eq.  (13)  of  r’]. 
Far-flelds 


For  the  far-fleld  conditions  given  by  Fqs.  (3.2-7)  and  (3.2-R), 

Eqs.  (3.8-2)  -  (3.8-4)  reduce  to  the  far-flelds  H=u^Hx  and  E=u1e..=u' r7M  where 

(f  tj)  a  fi '  a 

with  reference  to  Eq.  (4-62)  In  [25],  H.  Is  given  by 


j  1(0)  exp(-1kr)  P  cos(khcos9)  -  cos(kh) 


sln(kh)  27rr  I  sln9 


]  ,  n  £  q  < 

(3.8-q) 


7r/2 


0,  7T/2  <  9  <  7T 


The  time-averaged  radiated  power  density  s(fl)  =(^/2)|h^|  Is  given  by 


8(0)  = 


where 


^|l(0)|  f(0) 

87J^r^  s In^(kh) 


,  n  <  ft  <  7r/2 


(3.8-10) 


0,  7t/2  <  9  <  n 


f(e) 


cos( 


kh  cosft)  -  cos(kh) 


s  1  nft  j 

The  direction  of  maximum  radiated  power  Is  ft  =77/?  rad. 


The  Input  Impedance  of  a  monopole  element  of  length  h  on  a  groundplane  of 
Infinite  extent  Is  one-half  that  of  a  dipole  of  total  length  ?h  In  free  space 
Accordingly,  the  radiation  resistance  R  of  an  Infinitely  thin  element  on  a 
groundplane  of  Infinite  extent  Is  given  by  [compare  with  Fqs.  (4-70)  and 
(4-79)  In  [25]] 


D  -  ^  Cln(2kh)  + sln(2kh)[Sl(4kh)  -  2  Sl(2kh)l 

R  -  2  ^ 

47rsln  (kh) 

+  cos(2kh)[2  Cln(2kh)  -  Cln(4kh)]  ,  b=n 

^  (:i.R-ll) 


where  Sl(x)  and  Cln(x)  are  the  sine  Integral  and  modified  cosine  Integrals, 
respectively,  which  are  defined  following  Kq.  (3.2-32).  Kq.  (3.R-9)  Is  also 
approximately  valid  for  thin  elements  (kb  <<  1).  [For  example,  compare 
Eqs.  (3.2-33)  and  (3.2-35)  for  a  thin  monopole  element  In  the  absence  of  a 
groundplane. ] 

Substituting  Fqs.  (3.8-10)  and  (3.8-11)  Into  Fq.  (2.2-8),  the  numeric 
directive  gain  d(e)  Is  given  by 


A  f(e) 

B 


(3.8-1?) 


where 


B  =  Cln(2kh)  +  —  sln(2kh) [Sl(4kh)  -  2  Si(2kh)l 
+  cos(2kh)[2  Cln(2kh)  -  Cln(4kh)] 

Ff’r  a  quarterwave  monopole  element  (kh=“7r/2),  F,qs.  (3.8-11)  and  (3.8-1?) 


reduce  to 


47r 


Cln(270  ^  36.5395  ohms;  kh-7r/2,  b-0 


f(e) 


d(e)  = 


cos  [(7r/2)_  cosft]  ^ 


sln^G 


(3.R-13) 

(3.R-14) 


4  cos  [(7r/2)  cos^] 


sln^Q 


cin(2;r) 


,  kh=7r/2 


(3.«-15) 


The  peak  numeric  directivity  d(7r/2)  =  3.2818  corresponding  to  n(7r/7)  = 
10  d(7r/2)  =  5.1612  dBi.  The  directional  gain  is  plotted  In 

figure  8(f)  of  Section  3.9. 


For  an  electrically  short  monopole  element  (kh  «  1)  and 
)?  ss  120 TT  ohms ,  Eqs.  (3.8-11)  and  (3.8-12)  reduce  to 


R  =  10  (kh)^,  (kh)  «  1  (3.8-16) 

4 

f(e)  -  sln^e,  (kh)  «  1  (3.8-17) 

d(0)  3  sin  0,  (kh)  «  1  (3.8-18) 

The  peak  numeric  directivity  d(7r/2)  =  3.0  corresponding  to  rt(7r/2)  = 

10  log^Q  d(7r/2)  =  4.7712  dBi. 

Input  Impedance 

The  Input  resistance  R^^  Is  given  by  Eq.  (3.8-11)  which  is  exact 
for  an  Infinitely  thin  element  and  approximately  correct  for  thin 
elements  provided  that  the  element  current  distribution  Is  sinusoidal. 

The  input  reactance  X,  for  thin  elements  (kb  <<  1)  Is  given  by 

in 

[compare  with  Eqs.  (7-33)  and  (7-30)  in  [25) 
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Sl(2kh)  +  cos(2kh)[Sl(2kh)  - 


1 

2 


Si(4kh)] 


X 


In 


V 

47rsin^(kh) 


-  sln(2kh)[ln(h/b)  -  Cln(2kh)  +  Ctn(4kh)  +  Cln(kb^/h)] 

^  ^  (3.8-19) 


For  a  quarterwave  element,  In  the  limit  b— *0,  Eq.  (3.8-19)  gives  an 
input  reactance  =»  21.2576  ohms.  In  Eq.  (3.8-19),  for 

element  lengths  slightly  less  than  kh  =  (2n-l)(7r/2) ,  n»l ,  2,  3,  ..., 
which  are  approximately  one-half  the  resonant  lengths  for  a  monopole 
element  with  no  groundplane  [see  Eq.  (3.2-38)  ]. 


3.9  Summary  of  Results 


In  Section  3  a  sinusoidal  current  distribution  Is  assumed  on  the 
monopole  element.  Although  such  a  current  distribution  is  never 
exactly  realized  even  for  an  infinitely  thin  element  (see  Section  3.1), 
it  is  a  useful  approximation  for  sufficiently  thin  elements.  For 
example,  for  a  quarterwave  element  of  radius  b=10  ,  the  input 

Impedances  computed  by  determining  the  actual  element  current 
distribution  and  that  computed  by  assuming  a  sinusoidal  current 
distribution  differ  by  no  more  than  5%  for  groundplane  radii 
ka  «  6  -  14  (compare  table  7  of  Section  4.5  with  tables  A2-24  -  A2-42 
of  Appendix  A).  The  assumption  of  a  sinusoidal  element  current 
distribution  allows  for  models  which  are  computationally  simpler  in 
determining  input  Impedance  and  radiation  patterns  than  the  models 
which  follow  in  Section  4.  The  results  of  these  simpler  models  are 
summarized  in  this  section. 


The  electrical  properties  of  electrlcally-short  and  quarterwave 
monopole  elements  on  groundplanes  of  zero,  large,  and  infinite  extent 
are  compared  in  table  4  with  the  corresponding  properties  of 


Table  4 .  Electrical  Properties  of  Very  Thin  Monopole  Elements  on  Groundplanes  of  Zero 
Large,  and  Infinite  Extent 


electrlcally-short  and  half-wave  dipoles.  The  peak  directivity  of  a 
quarterwave  monopole  Is  1.88  dBl  and  5.16  dBl  for  groundplanes  of  zero 
and  infinite  extent,  respectively.  The  directivity  on  the  horizon  of 
a  quarterwave  monopole  is  1.88  dBl  and  -0.86  dBi  for  groundplanes  of 
zero  and  very  large  but  finite  extent,  respectively.  Slightly  smaller 
directivities  are  obtained  for  electrlcally-short  elements  than  for 
quarterwave  elements . 

The  radiation  resistances  obtained  by  different  methods  are 
compared  in  table  5  for  a  thin  quarterwave  element  on  a  groundplane  of 
radius  0  £  ka  £  8.5  wavenumbers.  The  suspected  best  available  results 
are  obtained  by  the  integral  equation  method  for  0  jC  ka  £  1.0,  by  the 
method  of  moments  (N»l)  for  1  <  ka  <  3.0,  by  the  oblate  spheroidal 
wave  function  method  for  3.00  ^  ka  ^  6.5,  and  by  the  method  of  moments 
(N“l)  for  6.5  <  ka  14.  The  results  obtained  by  the  method  of 
moments  is  in  good  agreement  with  the  suspected  best  available  results 
obtained  by  other  methods. 

The  numeric  directive  gain  patterns  of  a  quarterwave  element  on 
groundplanes  of  radii  ka  ■  0,  3,  4,  5,  yjhl ,  and  oo  wavenumbers  are 
plotted  in  figure  8.  These  polar  graphs  of  directive  gain  should  not 
be  confused  with  the  polar  graphs  of  radiated  power  density  plotted  in 
Ref.  [5]  for  constant  base  current.  It  should  be  noted  that  the  peak 
directivity  and  direction  of  peak  directivity  are  not  monotonlc 
functions  of  the  groundplane  radius. 

The  input  Impedance  of  a  thin  quarterwave  element  is  plotted  in 
figures  9  and  10  for  groundplane  radii  0  £  ka  £  14  wavenumbers.  The 
input  resistance  varies  between  19.4  and  46.1  ohms  and  asymptotically 
approaches  36.5  ohms  for  increasingly  large  groundplane  radii.  The 
input  reactance  varies  between  -  oo  and  +32.5  ohms  and  asjmiptotlcally 
approaches  +  21.3  ohms  for  increasingly  large  groundplane  radii. 
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Radiation  Resistance  (ohms) 


Normalised 
Groundplane  Radius  , 
2  It  a/A 

Integral  Equation 
Method 

Oblate  Spheroidal 
Wave  Function  Method 

Moment 

Method 

0 

*19.43 

0.25 

*19.48 

19.49 

0.50 

*19.62 

19.62 

0.75 

*19.86 

19.86 

1.00 

*20.23 

20.21 

1.25 

20.76 

*20.71 

1.50 

21.51 

*21.25 

1.75 

22.59 

*22.44 

2.00 

24.15 

*23.89 

2.25 

26.46 

*25.99 

2.50 

29.95 

27.32 

*29.02 

2.75 

35.44 

31.92 

*33.24 

3.00 

44.60 

*37.48 

38.62 

3.25 

*43.01 

44.12 

3.50 

*46.06 

47.57 

3.75 

*45.55 

47.35 

4.00 

*42.67 

44.43 

4.25 

*39.23 

40.58 

4.50 

*36.23 

37.13 

4.75 

*34.00 

34.46 

5.00 

*32.57 

32.68 

5.25 

*31.93 

31.70 

5.50 

*32.13 

31.53 

5.75 

*33 . 23 

32.26 

6.00 

*35.23 

34.04 

6.25 

*37.85 

36.94 

6.50 

*40.33 

40.56 

6.75 

30.12 

*43 . 53 

7.00 

30.09 

*44.20 

7.25 

*42.30 

7.50 

*39.10 

7.75 

*35.96 

8.00 

*33.50 

8.25 

*31.88 

8.50 

*31.16 

*SDSPECTED  BEST  AVAILABLE  RESULT 


Table  5.  Radiation  Resistance  of  a  Thin  Quarterwave  Eleirent  at 
the  Center  of  a  Circular  Groundplane  of  Radius 
0<ka<8.5  Wavenumbers.  (Sinusoidal  current  distribution 
assumed  on  element.) 
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Input  Reactance  of  a 
Groundplane  of  Radius 


r/d 


The  directive  gain  on  the  horizon,  peak  directivity,  and  elevation 
angle  of  the  peak  directivity  of  a  quarterwave  element  are  plotted  in 
figures  11  -  13  for  groundplane  radii  0  £  ka  £  14  wavenumbers.  The 
directive  gain  on  the  horizon  varies  between  +1.9  and  -1.9  dBl  and 
asymptotically  approaches  -0.86  dBl  for  increasingly  large  but  finite 
groundplane  radii.  The  peak  directivity  varies  between  +1.0  and 
+5.3  dBl  and  asymptotically  approaches  +5.2  dBi  for  increasingly  large 
groundplane  radii.  The  elevation  angle  of  the  peak  directivity  varies 
between  33  and  90  degrees  and  asymptotically  approaches  90  degrees  for 
increasingly  large  groundplane  radii. 

The  radiation  resistances  of  thin  elements  of  length  h/\»0.25, 

0.1  and  0.025  for  groundplane  radii  0  £  ka  £  8  wavenumbers  are 
compared  in  figure  14.  The  radiation  resistance  of  each  element  is 
normalized  to  its  radiation  resistance  for  a  groundplane  of  zero 
extent.  The  normalized  radiation  resistance,  as  a  function  of 
groundplane  radius,  is  approximately  independent  of  the  element 
length.  For  an  electrically-short,  thin  element  (h/X  whose 

length  is  small  or  comparable  to  the  groundplane  radius,  the  input 
reactance  (not  shown)  is  approximately  independent  of  groundplane 
radius. 


Directive  Gain  on  Horizon  of  Thin  Quarterwave  Element  at  the  Center  of  a  Circular 
Groundplane  of  Radius  0<k.a<14  Wavenumbers 


Figure  12  .  Peak  Directivity  of  a  Thin  Quarterwave  Element  at  the  Center  of  a  Circular 
Groundpla.nB  of  Rndius  O^k.3.^14  WavenuinbBirs 


Elevation  Angle  of  Peak  Directivity  for  a  Thin  Quarterwave  Element  at  the  Center 
of  a  Circular  Groundplane  of  Radius  0^ka<14  Wavenumbers 


liation  Resistance  of  Various  Length  Thin  Elements  at 
rcular  Groundplane  in  Free  Space 


SECTION  4 


MODELS  IN  WHICH  THE  CURRENT  DISTRIBUTIONS  ON  THE  MONOOOLE 
ELEMENT  AND  GROUNDPLANE  ARE  ROTH  INITIALLY  UNKNOWN 


4.1  Boundary  Conditions 

In  Section  3  the  total  field  at  an  arbitrary  field  point  could  be 
expressed  simply  as  the  vector  sum  of  an  Incident  field  and  an  Induced 
field  [see  Eq.  (3.1-12)1  because  the  Incident  field  was  specified  from 
the  Initially  known  current  distribution  on  the  monopole  element.  In 
this  section,  such  a  procedure  Is  not  possible  because  the  current 
distributions  on  the  element  and  groundplane  are  both  Initially 
unknown.  Instead,  In  this  section,  the  total  field  Is  determined  by 
representing  the  unknown  current  distribution  on  either  the  element, 
groundplane,  or  both  by  a  finite  series  of  overlapping  modes  with 
unknown  current  amplitudes.  The  current  amplitudes  are  determined  by 
matrix  Inversion  subject  to  boundary  conditions  comprising  current 
constraints  and  the  excitation  voltage  across  the  coaxial  line  feed. 

The  antenna  geometry  Is  shown  In  figure  I  of  Section  2.1.  The 
current  waveform  Is  given  by  Eq.  (2.1-2).  The  models  which  follow  In 
this  section  are  based  on  the  current  characterization  and  circuit 
Idealization  In  figure  2(c)  of  Section  2.4. 

The  coaxial  line  feed  and  excitation  voltage,  which  Is  not 
explicitly  shown  In  figure  1,  Is  characterized  In  figure  2(c)  by  a 
surface  magnetic  current  density  (magnetic  frill)  M^  which  sits  on  top 
of  a  thick  groundplane  of  radius  b^  where  b^  Is  the  radius  of  the 
outer  conductor  of  the  coaxial  line  feed.  The  magnetic  frill  Is 
defined  over  the  groundplane  by  Fq.  (2.4-R). 


ss 


Constraints  on  the  various  circuit  currents  of  figure  2(c)  are 
given  by  Eqs .  (2.4-3)  -  (2.4-7).  These  constraints,  together  with  the 
magnetic  frill,  constitute  the  boundary  conditions  on  the  current 
amplitudes . 

4.2  Method  of  Moments,  0  <  ka  <14 

When  the  element  and  groundplane  current  distribution  are 
initially  unknown,  the  current  distributions  may  be  determined  by  a 
sinus oldal-Galerkln  moment  method  employed  by  Richmond^^\  The 
antenna  geometry  is  shown  in  figure  1  of  Section  2. 

In  the  moment  method,  the  element  is  divided  into  N  equal  segments 
(see  figure  15)  of  length  d'  given  by 

d'  ”  h/N,  N  is  a  positive  integer  (4.2-1) 

where  h  is  the  element  length.  The  groundplane  is  divided  into  M 
concentric  annular  zones  (see  figure  16)  of  width  d  given  by 

d  »  (a-b)/M  >  bj^  -  b,  M  is  a  positive  integer  (4.2-2) 

where  a  is  the  groundplane  radius ,  b  is  the  element  radius ,  and  b^  is 
the  radius  of  the  outer  conductor  of  the  coaxial  line  feed.  The 
groundplane  extends  from  the  coaxial  line  inner  conductor  because  of 
the  equivalent  circuit  representation  in  figure  2(b)  of  the  coaxial 
line  excitation.  The  element  current  distribution  I(z)  and 
groundplane  current  distribution  I(P)  are  the  sum  of  the  current 
distributions  on  each  element  segment  and  groundplane  annular  zone, 
respectively,  and  are  given  by 


Figure  15.  Element  with  N  Equal  Segments  of  Length  d' 
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(4.2-4) 
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K")  •  L  i;(p) 

m“l 

where  Ijj(z)  and  I'(P)  are  the  current  distributions  on  the  n*"^  segment 
and  annular  zone,  respectively. 


In  the  sinusoidal  Galerkln  method,  the  current  distribution  I  (z) 

th  ^ 

on  the  Ti  element  segment  Is  approximated  by 


In(z)  « 


sln[k(nd'-z)]  +  8ln{k[z-(n-l)]d'} 


sin(kd') 


(n-l)d'  ^  z  £  nd';  n*l,2,3,  ....N 

(4.2-5) 


where  *  l(z=0)  and  I 


N+1 


0. 


The  current  distribution  1^(P)  on  the  annular  zone  is  approximated 
by 


K(P) « 


m 


r  8ln{k[md+b-p]}  +  8ln{ktp-(m-l)]d-b} 

sin(kd)  * 

b  +  (m-1  )d  £  P  £  b+md;  m>*l ,  2 ,  3,  . .  .M 

(4.2-6) 


where  I'  .  *  -I  .  ■  -l(z*0)  and  I 
m-1  n-1  '  ' 


M+1 


The  current  distribution  on  each  segment  and  annular  zone  is 
therefore  the  sura  of  two  overlapping  dipole  modes  except  for  the 
segment  and  annular  zone  adjacent  to  the  base  of  the  element.  The 
"base"  mode  may  be  thought  of  as  a  dipole  mode  comprising  an  element 
monopole  mode  (n=l)  and  a  groundplane  monopole  mode  (m=l)  with 
terminals  along  the  circumference  of  the  base  of  the  element.  The 
element  and  groundplane  dipole  modes  are  functionally  tied  to  each 
other  by  means  of  the  base  mode  amplitude  constraint  " 

I(z-O). 
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Hie  element  and  groundplane  total  current  distributions  are 

represented  In  Eqs .  (A. 2-3)  -  (A. 2-6),  as  a  series  of  N+M  overlapping 

sinusoidal  dipole  modes  with  W-M  unknown  complex  current  amplitudes. 

The  N+M  amplitudes  are  determined  In  Richmond's  moment  method  by 

Inversion  of  a  (N+M)  x  (N+M)  matrix  subject  to  the  boundary  conditions 

discussed  in  Section  A.l.  The  constraint  I  ,=  -I'  ,  reduces  the 

n=»l  m^l 

matrix  size  to  (N4M-1)  x  (N+M-1). 

The  numbers  of  subdivisions,  N  and  M,  are  limited  by  the  cost  of 
computation  time  and  by  the  precision  of  the  computer.  The  accuracy 
of  the  solution  can  be  decreased  appreciably  if  either  N  or  M  is  too 
small  or  too  large.  The  method  of  moments  converges  to  a  solution 
when  an  increase  or  decrease  of  unity  in  the  value  of  N  or  M  does  not 
appreciably  alter  the  solution  for  input  impedance. 

A  method  of  moments  computer  program  for  a  monopole  element  at  the 
center  of  a  circular  groundplane  in  free  space  was  obtained  from 
Prof.  Richmond  of  Ohio  State  University.  The  program  computes  the 
input  Impedance  and  the  N+M  complex  current  amplitudes  on  the  element 
and  groundplane  for  a  voltage  input  V(0)  1  volt.  The  program, 

written  in  FORTRAN  IV  language  and  in  single  precision,  was  edited  and 

converted  by  The  MITRE  Corporation  to  double  precision  for  use  on  a  DEC 

VAX  11/780  computer.  A  listing  of  the  MITRE  version,  designated 
"RICHMDl",  is  given  in  Appendix  B2. 

Our  experience  with  this  program  suggests  the  following 
constraints  on  the  use  of  the  program.  At  least  double  precision  is 
required  to  give  convergent  results.  Meaningful  results  were  obtained 
for  element  radii  b/ X  >10  and  groundplane  radii  ka  >  0.25.  As  a 
rule  of  thumb,  N*2-3  kh  and  M=2-3  ka  give  reasonably  accurate  and 
convergent  results.  The  amount  of  central  processing  unit  (CPU)  time 
on  the  VAX  11/780  computer  is  approximately  2  minutes  for  N+M  =  20,  6 


E 

I 


f 


I 


i 


minutes  for  N+M=30,  34  minutes  for  N+-M=60,  and  50  minutes  for  N+M“75. 
It  Is  most  likely  that  the  CPU  time  could  be  considerably  reduced  If  a 
more  efficient  method  than  the  Grout  method  were  employed  for  matrix 
inversion. 

The  versatility  of  a  sinusoidal  -  Galerkln  method  of  moments  Is 
illustrated  by  the  case  N*!.  A  sinusoidal  current  distribution  Is 
imposed  on  the  element  by  setting  N=1  [see  Eq.  (4.2-5)]. 

The  element  and  groundplane  current  distributions  obtained  from 
program  RICHMDl  can  be  utilized  to  obtain  the  far-fleld  radiation 
pattern  In  the  following  manner. 

The  magnetic  field  Intensity  H(P)  at  a  far-fleld  point  P(r,  0,  0) 
and  for  an  e^*^  waveform  Is  related  to  the  magnetic  vector  potential 
T(P)  by^^^^ 


H(P)  “  -j(l/^)WU^  X  A(P) 


(4.2-7) 


where  u^^  is  a  unit  vector  in  the  radial  direction  and  ^  is  the 

free-space  wave  Impedance.  The  magnetic  vector  potential  A(P)  at  the 

far-fleld  point  P  resulting  from  a  current  source  point  Q(x',  y',  z') 

=•  Q(b,  b,  z')  on  the  element  or  a  source  point  Q(P',  7T/2,  0')  on  the 

(  26  ) 

groundplane  (see  figure  4)  is  given  by'' 

MP)  -  (/^Q/47r)  J J(Q)  exp(-jk^)  (1/^)  dv 

V 

«  [jU  exp(-jkr)/47Tr]  f  J(Q)  exp(jku^-OQ)  dv  (4.2-8) 


where 


PQ  -  tr^-2r'u^*0Q  +  x'^+y'^+z'^] asT-  ^  »  (x'^+y'^+z'^)^^^ 

u  •  00  =  x'  sin0  cos^  +  y'  sin0  sin(ll  +  z'  cos0 
r  _7 

=  permeability  of  free  space  =  47r  x  10  henry/m. 

J(Q)  =  source  surface  current  density  at  an  arbitrary  source  point  Q 

on  either  the  element  or  groundplane  (amp/m) 

r  =  OP  *  radial  distance  between  the  origin  0(0, 0,0)  and  the 

far-field  point  P(r,0  ,0) 

dv  =  differential  surface  area  containing  the  current  source 
points  Q  (m2) 

Substituting  Eq.  (4.2-3)  into  Eq.  (4.2-8),  the  far-field  magnetic 

C  6  ) 

field  Intensity  resulting  from  the  monopole  current  distribution 

la  given  by 


(e)  j  exp(-jkr)  sln0 


J  exp(-jkr) 

47n:  8in0  8ln(kd') 


■/  J,'- 


r  ” 


(z')  exp( jkz'co80)  dz' 


exp( jnkd'cos0) 


^l-cos(kd')  cos(kd'co30)  -  cos  0  sln(kd' )8ln(kd'cos0) 

+  J  cos(kd')  sin(kd'  cos0)  -  j  cos  0  8in(kd')  cos(kd'cos0)j 

exp[j(n-2)  kd'cos0]J  |^co8(kd' )  cos(kd'co8  0) 

+  cos  0  sln(kd')  sln(kd'cos  0)  -  1  +  j  cos(kd')  sln(kd'cos  0) 
-  J  cos  0  sln(kd')  co3(kd'cos  0)j  1  (4.2-9) 


•V  * 


Substituting  Eq.  (4.2-4)  Into  Eq.  (4.2-8),  the  far  field  magnetic 

(a) 

field  Intensity  Hjjj  resulting  from  the  groundplane  current 
distribution  Is  given  by 


^  a  -  27r 

•'o  •'0  “=1 


cos  0  cos((ll'-(|))  ddl'  dP' 


(4.2-10) 


Eq.  (4.2-10)  Is  not  readily  evaluated  when  Eq .  (4.2-6)  Is  substituted 
Into  it.  An  approximate  simplified  solution  can  be  found  when  each 
annular  zone  of  the  groundplane  Is  subdivided  Into  X  smaller  annular 
zones  of  width  AP  =  d/X  so  that  the  current  distribution  In  each 
subdivided  zone  Is  almost  constant  and  approximately  equal  to  Its 
average  current.  The  total  number  J  of  subdivided  annular  zones  on 
the  entire  groundplane  Is  given  by 


J  =  XM 


(4.2-11) 


The  average  current  In  the  uth  subdivided  zone  Is  given  by 


c  f 


u  APt-b 


sln[k(md+b-p]  +  I' ^  ^  sln{k[p-  (m-l)d  -  b] } 


)AP+  b 


sln(kd) 

(W-(m-l)d]<P<(b+ind) 


(4.2-12) 


The  current  distribution  I(P)  on  the  groundplane  Is  therefore  given 
approximately  by 


I(P)  «  yi  I 


(4.2-13) 


Substituting  l^(P')  for  I^(P')  in  Eq.  (4.2-10), 


a  27r 


“(j)  « 


-j  exp(-jkr) 
Attt 


// 


^  I^(P')exp(  jk,/:j'cos((>'-0)  ] 


0  0 


^1=1 


cos9  cos((|l'-(|))d()'dp' 


^xp(- jk.r)cos  Q 
Attt  sin0 


IJP')  fjoik  sine[(u-l)AP +b]}  -  J^{ksine[u  A/^+b]  }1 
^  (4.2-14) 


where  J  is  the  Bessel  function  of  the  first  kind,  of  order  zero, 
o 


The  total  far-field  magnetic  field  Intensity 


is  given  by 
(4.2-15) 


fe)  Ce) 

where  "^and  H^®^are  given  by  Eqs .  (4.2-9)  and  (4.2-14),  respectively. 


The  time-averaged  radiated  power  density  s(0),  directive  gain 
d(0),  and  radiation  resistance  R  are  found  from  Eqs.  (2.2-4), 


and  (2.2-7)  where  H 


(total) 


is 


given  by  Eq.  (4.2-15)  and  l(z='0)  =  ^n=l' 


Numerical  evaluation  of  s(0),  d(0),  and  R  are  implemented  by  means  of 
MITRE  computer  program  "RICE1MD2"  written  in  FORTRAN  language  for  use 
on  a  DEC  VAX  11/780  computer.  A  listing  of  program  RICHMD2  is  given 
In  Appendix  B2. 


A  discussion  of  the  results  obtained  by  the  method  of  moments,  for 
a  sinusoidal  current  distribution  on  the  element  (N=l),  is  given  in 
Section  3.4. 


Computer  printouts  of  the  directive  gain  and  far-field  radiation 
pattern  of  the  experimental  monopole  antennas  discussed  in  Section  5 
were  obtained  by  the  method  of  moments  and  are  given  in  Appendix  A2. 


The  input  Impedance  determined  by  the  method  of  moments  (N=4, 
M»3ka)  for  a  quarterwave  element  of  radius  b“10  and  groundplane 
radii  6  ^  ka  30  wavenumbers  is  compared  in  table  7  of  Section  4.5 
with  that  obtained  by  the  method  of  moments  combined  with  the 
geometric  theory  of  diffraction  (GTD).  For  ka  >  14,  Richmond's  method 
of  moments  is  Inaccurate.  In  addition  to  the  large  CPU  time  required 
(greater  than  30  minutes  for  ka  >  14),  convergent  results  are 
difficult  to  obtain  for  ka  >  14.  It  is  most  likely  that  the 
usefulness  of  Richmond's  method  of  moments  could  be  extended  to 
significantly  larger  groundplane  radii  if  a  more  efficient  method  of 
matrix  inversion,  different  from  the  Crout  method  employed  by 
Richmond,  were  utilized.  When  the  element  current  distribution  is  not 
assumed,  Richmond's  method  of  moments  presently  gives  the  best 
available  results  for  groundplane  radii  ka  <  14  wavenumbers. 

The  lower  limit  of  ka,  for  accurate  results  utilizing  Richmond's 
method  of  moments,  has  not  been  firmly  established.  For  an  element 
segment  number  N  >  1,  we  have  obtained  useful  results  for  ka  as  small 
as  0.25  wavenumbers  but  not  for  groundplanes  of  zero  extent.  For  N=*l, 
accurate  results  were  obtained  (subject  to  the  constraint  of  a 
sinusoidal  element  current  distribution)  for  groundplanes  of  zero 
extent . 

The  radiation  resistance  of  various  diameter  resonant  (zero  input 
reactance)  elements  is  compared  in  figure  17  of  Section  4.5  with  the 
radiation  resistance  of  an  Infinitely  thin  quarterwave  element  for 
groundplane  radii  0  £  ka  _<  8.5  wavenumbers.  For  ka  — ►O,  the 
radiation  resistance  of  the  resonant  elements  is  appreciably 
different  from  that  of  the  infinitely  thin  quarterwave  element  because 
an  infinitely  thin  element  on  a  groundplane  of  zero  extent  is  resonant 
for  an  element  length  equal  to  0.5X  and  not  0.25 X  [see  Eq .  (32-38)]. 

The  lengths  and  radiation  resistances  of  various  diameter  resonant 
elements,  for  groundplanes  of  radii  0.25  <  ka  7  wavenumbers,  are 
plotted  in  figures  18  and  19  of  Section  4.5. 


A. 3  Method  of  Moments  Combined  with  Geometric  Theory  of  Diffraction, 
8.5  <  ka  <00 


The  method  of  moments ,  when  used  to  determine  the  current 
distributions  on  both  the  element  and  the  groundplane,  can  require 
considerable  computation  time  for  large  groundplane  radii.  For 
example,  Richmond's  method  of  moments  requires  approximately  one-half 
hour  of  CPU  time  on  the  VAX  11/780  computer  for  a  groundplane  radius 
ka=20  with  N=4  element  segments  and  M=60  groundplane  annular  zones 
(see  Section  4.2).  Although  it  may  be  possible  to  reduce  the 
computation  time  by  use  of  a  more  efficient  program  than  the  one 
employed  by  Richmond,  it  is  of  Interest  to  find  an  alternative  method 
for  large  groundplanes . 

Awadalla  and  Maclean^^ ^ ^ have  reduced  the  computation  time  for 
large  groundplanes  by  combining  the  method  of  moments  with  the 
geometric  theory  of  diffraction  (GTD).  The  element  current 
distribution  is  determined  by  the  method  of  moments  and  the  effect  on 
input  Impedance  by  groundplane  edge  diffraction  is  determined  by  GTD. 

The  antenna  geometry  is  shown  in  figure  1  of  Section  2.1.  The 
method  of  Awadalla  and  Maclean  is  described  in  [9]  for  determining  the 
input  impedance  and  in  [10]  for  determining  the  radiation  pattern. 

A  fictitious  magnetic  edge  current  is  defined  and  expressed  in 
terms  of  a  GTD  diffraction  coefficient  to  account  for  diffraction  by 
the  edge  of  the  groundplane.  The  method  of  GTD  is  valid  only  for 
sufficiently  large  groundplane  radii.  Unfortunately,  this  method,  as 
applied  by  Awadalla  and  Maclean,  does  not  determine  the  groundplane 
current  distribution. 
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The  groundplane  current  distribution  Is  Idealized  by  Avadalla  and 
Maclean  in  determining  the  radiation  pattern.  The  current 
distribution  on  the  top  of  the  groundplane  is  assumed  to  be  that  for 
an  infinite  groundplane  but  defined  only  over  the  finite  extent  of  the 
groundplane.  The  net  current  at  the  edge  of  the  groundplane  is  set 
equal  to  zero.  This  constraint  satisfies  the  boundary  condition  given 
by  Eq.  (2.4-4)  and  requires,  at  the  edge  of  the  groundplane,  that  the 
current  on  the  bottom  of  the  groundplane  be  equal  but  in  the  opposite 
direction  to  that  on  the  top  of  the  groundplane.  The  current 
distribution  on  the  bottom  of  the  groundplane  is  then  assumed  to 
decrease  linearly  from  the  edge  of  the  groundplane  to  zero  at  the  base 
of  the  element . 

The  method  of  moments  employed  by  Awadalla  and  Maclean  is  not  a 
sinusoldal-Galerkin  method.  Consequently,  one  cannot  impose  a 
sinusoidal  current  distribution  on  the  element  by  setting  the  number 
of  element  segments  N-1  as  a  test  case  for  purposes  of  comparison  with 
other  models. 

A  program  listing  was  obtained  from  Prof.  Awadalla  and  edited  at  The 
MITRE  Corporation.  A  listing  of  the  edited  program  "AWADALLA”  is 
given  in  Appendix  A5.  Program  AWADALLA  is  written  in  FORTRAN  language 
for  use  on  a  DEC  VAX  11/780  computer.  The  CPU  time  for  N*»30  element 
segments  is  less  than  10  seconds. 

Program  AWADALLA  was  utilized  to  obtain  the  input  impedance, 
directive  gain,  and  far-fleld  elevation  pattern  of  a  quarterwave 
element  of  radius  b=10  on  groundplane  of  radii  8  £  ka  £  50 

wavenumbers.  The  results  for  input  Impedance  are  given  In  table  6. 

These  results  are  compared  in  table  7  of  Section  4.5  with  those 
obtained  by  the  method  of  moments  (N=»4).  Concerning  input  Impedance, 
the  method  of  moments  combined  with  GTD  gives  useful  results  for 
ka  >  6,  accurate  results  for  ka  >  8,  and  the  suspected  best  available 
results  for  14  <  ka  < oo  . 


TABLE  6.  INPUT  IMPEDANCE  OF  A  QUARTERWAVE  ELEMENT  OF  RADIUS  ly-10"^X  AT  THE 
CENTER  OF  A  CIRCULAR  GROUNDPLANE  OF  RADIUS  8  <  ka  <  50  WAVENUMBERS 

(Method  of  moments  combined  with  geometric  theory  of  diffraction) 


INPUT  RESISTANCE 

INPUT  REACTANCE 

a 

(OHMS) 

(OHMS) 

A  sample  computer  printout  of  directive  gain  and  the  far-field 
elevation  pattern  obtained  by  program  AWADALLA  for  ka“20  is  given  in 
Appendix  A5.  Unfortunately,  the  Idealizations,  made  by  Awadalla  and 
Maclean  in  characterizing  the  groundplane  current  distribution,  yield 
unrealistic  peak  directivities  and  elevation  patterns.  For  example, 
for  ka»49,  a  peak  directivity  of  7.5  dBl  was  obtained.  This  result 
for  a  thin  quarterwave  element  seems  unlikely  because  a  thin 
quarterwave  element  with  a  sinusoidal  current  distribution  and  mounted 
on  a  groundplane  of  infinite  extent  has  a  peak  directivity  of  only 
5.2  dBl. 

The  method  of  moments  combined  with  GTD  via  the  definition  of  a 
fictitious  magnetic  edge  current  has  also  been  reported  by  Thiele  and 
Newhouse^^^^  for  computing  input  impedance  and  by  Stutzman  and 

Thlele^^®^  for  computing  the  far-field  radiation  pattern.  In  the 

f  28^ 

method  of  Stutzman  and  Thiele'  ' ,  the  far-field  pattern  is  determined 
without  idealizing  the  groundplane  current  distribution. 

4. A  Method  of  Images,  ka=*  oo 

For  the  idealized  case  of  a  monopole  element  mounted  on  a 
groundplane  of  infinite  extent  and  of  infinite  conductivity,  a 
monopole  antenna  of  length  h  may  be  modelled  by  the  method  of  images 
as  a  dipole  of  total  length  2h  but  with  one-half  the  input  impedance 
and  double  the  peak  numeric  directivity  of  the  dipole  (see 
Section  3.8).  The  infinite  groundplane  prevents  radiation  into  the 
hemisphere  below  the  groundplane  but  generates  fields  in  the  upper 
hemisphere  identical  to  those  of  a  dipole. 

A  detailed  treatment  of  the  fields  and  input  impedance  of  a  dipole 
is  given  in  [13].  An  excellent  summary  of  the  present 
state  of  the  art  of  dipole  theory.  Including  plots  of  input  Impedance 
as  a  function  of  dipole  length  and  dipole  radius,  is  given  in  [14]. 


4.5  Summary  of  Results 


In  Section  4,  the  current  distribution  on  the  element  is 

determined  by  the  method  of  moments  rather  than  being  assumed  as  was 

the  case  in  Section  3.  The  determination  of  the  element  current 

-4  ^ 

distribution  is  particularly  important  for  element  radii  b  >  10  A  . 

The  essential  difference,  in  the  two  models  utilized  in 
Section  4,  is  the  treatment  of  the  groundplane  current  distribution. 

In  Richmond's  model,  the  groundplane  current  distribution  is 
determined  by  the  method  of  moments.  Useful  results  are  obtained  for 
0  <  ka  £  14  wavenumbers.  In  the  model  of  Awadalla  and  Maclean  for 
large  but  finite  groundplanes ,  the  input  Impedance  is  accurately 
computed  by  the  introduction  of  a  fictitious  magnetic  edge  current 
determined  by  the  Geometric  Theory  of  Diffraction  (GTD).  However,  in 
that  model  the  groundplane  current  distribution  is  not  determined  but 
Instead  is  idealized  —  causing  unrealistic  peak  directivities  and 
far-field  radiation  patterns. 

The  input  impedances,  determined  by  the  method  of  moments  and  the 
method  of  moments  combined  with  GTD,  are  compared  in  table  7  for  a 
quarterwave  element  of  radius  b»10  ^A  and  groundplane  radii 
6  £  ka  £  30  wavenumbers.  The  suspected  best  available  results  for 
input  impedance  are  obtained  by  the  method  of  moments  for  0  <  ka  £  8.5 
and  by  the  method  of  moments  combined  with  GTD  for  8.5  <  ka  <  oo  .  The 
method  of  moments  combined  with  GTD  is  Inaccurate  for  ka  ^  8 
(approximately  6Z  error  in  input  reactance  for  ka“8  and  7Z  error  in 
input  resistance  for  ka-6).  Richmond's  method  of  moments  is 
inaccurate  for  ka  >  14  (approximately  IIZ  error  in  input  reactance  and 
5%  error  in  input  resistance  for  ka=15).  Richmond's  method  of  moments 
would  be  useful  for  ka  >  14  wavenumbers  if  a  more  efficient  method 
than  the  Grout  method  were  employed  for  matrix  inversion. 
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The  radiation  resistance  of  various  diameter  resonant  (zero  Input 
reactance)  elements  Is  compared  In  figure  17  with  the  radiation 
resistance  of  an  Infinitely  thin  quarterwave  element  for  groundplane 
radii  0  ^  ka  £  8.5  wavenumbers.  The  radiation  resistance  of  the 
resonant  elements  is  not  too  appreciably  different  from  the 
quarterwave  element  for  ka  >  1  wavenumber.  However  as  ka  — ►O,  the 
radiation  resistances  of  the  resonant  elements  become  increasingly 
different  from  that  of  the  infinitely  thin  element  because  an 
infinitely  thin  element  on  a  groundplane  of  zero  extent  Is  resonant 
for  an  element  length  equal  to  0.5X  and  not  0.25X  [see  Eq.  (3.2-38)]. 

The  lengths  and  radiation  resistance  of  various  diameter  resonant 

elements  for  groundplane  of  radii  0.25  ^  ka  ^  7  wavenumbers  are 

plotted  in  figures  18  and  19.  For  these  groundplane  radii,  the 

element  resonant  length  h  / X  varies  from  approximately  0.22  to  0.34 

—7  —2 

wavelengths  for  element  radii  10  ^  b/X  <  10  wavelengths.  The 

resonant  radiation  resistance,  for  these  groundplane  and  element 
radii,  varies  from  approximately  21  to  65  ohms. 
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Figure  18  .  Resonant  Length  of  Various  Diameter  Elements  at  the  Center  of  a  Circular 
Groundplane  of  Radius  0.25<ka<7  Wavenumbers 


SECTION  5 


COMPARISON  WITH  EXPERIMENTAL  RESULTS 


The  Input  Impedance  and  elevation  patterns  of  several  monopole 
antennas  were  measured,  each  at  a  different  frequency  within  the 
frequency  band  30-25A  MHz,  on  the  MITRE  Corporation  VHP  antenna  range. 
Each  antenna  consisted  of  a  tubular  element  of  radius  b=0.25  In. 
mounted  at  the  center  of  a  groundplane  of  radius  a=4  ft.  and  fed  by  a 
50  ohm  RG-214  coaxial  cable  with  a  type  N  panel  connector  and  a  50  ohm 
tapered  adapter  to  the  element.  The  length  h  of  each  element  was 
approximately  a  quarterwave.  The  exact  length  of  each  element  length 
was  chosen  to  be  resonant  (Input  reactance  =  zero  ohms)  for  a 
groundplane  of  infinite  extent.  On  a  groundplane  of  finite  extent, 
the  Input  reactance  is  expected  to  asymptotically  approach  zero  as 
ka  — ►  00  . 

The  VHP  antenna  range  Is  located  on  the  roof  of  MITRE  E-Buildlng. 
The  transmitting  and  receiving  antennas  are  at  a  height  27  ft.  above 
the  roof  and  are  separated  by  approximately  40  ft.  A  conducting  fence 
16  ft.  high  and  48  ft.  wide  is  located  on  the  roof  midway  between  the 
transmitting  and  receiving  antennas  to  minimize  multipath  reflections 
from  the  roof.  Lossy  ferrite  toroids  (Perronlcs,  Type  B  material, 
0.540  In.  I.  D.  X  0.870  in  0.  D.  x  0.25  in.)  are  spaced  along  the 
receiving  and  transmitting  cables  to  minimize  currents  on  the  outside 
of  the  cables.  Outside  rf  Interference  Is  reduced  to  at  least  40  dB 
below  the  desired  signal  by  the  use  of  narrowband  rf  filters. 


The  monopole  test  antenna  was  operated  In  the  receiving  mode  to 
obtain  elevation  patterns  and  on  the  same  platform  for  Input  Impedance 
measurements.  The  monopole  groundplane  was  supported  by  a  10  ft. 
wooden  vertical  mast  mounted  on  an  antenna  pedestal. 

The  input  Impedance  was  measured  with  a  Hewlett-Packard  HP-8754A 
Network  Analyzer  with  a  computerized  printout.  The  measurement  test 
set-up  is  shown  in  figure  20.  The  predicted  input  Impedance  of  each 
experimental  monopole  antenna  was  determined  by  Richmond's  method  of 
moments  utilizing  program  RICHMDl.  The  predicted  and  measured  input 
Impedances  are  compared  in  Table  8.  The  measured  input  resistance 
differs  from  predicted  values  by  approximately  1-10%  over  a  range  of 
normalized  groundplane  radii  =  0.77  -  6.5  wavenumbers.  The  measured 
input  reactance  differs  by  approximately  2-12  ohms  from  the  predicted 
values.  Since  some  of  the  predicted  values  are  near  resonance,  a 
percentage  comparison  la  not  made  for  input  reactance. 

Measured  elevation  patterns  of  most  of  the  test  monopole  antennas 
specified  in  Table  8  are  compared  in  figures  21  -  29  with  theoretical 
patterns  predicted  by  the  method  of  moments  (see  Appendix  A2  for  the 
computer  printouts).  Allowing  for  some  multipath  distortion  by  the 
VHF  test  range,  the  received  patterns  are  in  good  agreement  with  the 
predicted  pattern.  It  should  be  noted  the  measured  pattern  for  ka  = 
0.766  is  not  appreciably  different  from  that  predicted  for  the 
monopole  element  Itself.  The  effect  of  the  groundplane  on  the  shape 
of  the  pattern  is  not  readily  apparent  until  ka  >  2.0. 
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REFERENCE  PLANE  EXTENSION  CABLE 
OF  TWICE  TKE  LENGTH 
OF  TEST  CABLE 


HP2673A 


PRINTER 


BASIC  LANGUAGE 

BASIC  LANGUAGE  EXTENSIONS 


MITRE  PROGRAM  'NT WK ANAL' 


Impedance  of  0.50  Diameter  Monopole  Elements  on  an  8  ft.  Diameter 
plane.  (Comparison  of  moment  method  predictions  with  MITRE  measurements. 


Measured  and  Theoretical  Elevation  Patterns,  30  MHz 


retical  Elevation  Pat 


Measured  and  Theoretical  Elevation  Patterns,  43  MHz 


Measured  and  Theoretical  Elevation  Patterns,  54  MHz 


Measured  and  Theoretical  Elevation  Patterns 


Measured  and  Theoretical  Elevation  Patterns,  75  MHz 
(h/A  =0.2374,  b/A  =  15.89  x  10"'^,  ka  =  1.915) 


Measured  and  Theoretical  Elevation  Patterns,  86  MHz 


Measured  and  Theoretical  Elevation  Patterns,  117  MHz 


Measured  and  Theoretical  Elevation  Patterns,  156  MHz 
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